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DISC-O-TIC 

A  Discrete-Time  Analytical  Meta-Model 
for  Use  in  Combat  Systems  Studies 
that  Utilize  High-Resolution  Simulation  Models 


Patricia  A.  Jacobs 
Donald  P.  Gaver 


EXECUTIVE  SUMMARY 

This  paper  provides  various  models  for  extending  or  extrapolating  in  time,  and 
varying  and  enhancing  in  coverage  capability,  the  output  of  a  high-resolution  simulation 
model,  here  specifically  but  not  exclusively,  the  Army’s  COSAGE.  The  models  we 
propose,  generically  called  DISC-O-TIC,  which  is  short  for  “Discrete-Time  Analytical 
Meta-Model)  are  tailored  to  employ  the  discrete-time  output  of  COSAGE  and, 
potentially,  many  other  such  models.  The  state  variables  of  the  two  opponent  forces  are 
the  numbers  of  vehicles  or  platforms  of  a  variety  of  types  for  both  sides,  measured  at 
regular  (discrete)  times  only:  e.g.  12-hour  intervals,  1-day  intervals,  2-day  intervals,  etc. 
as  selected  by  the  analyst.  The  outputs  are  survivors  of  the  battle  and  munitions  expended 
after  a  given  number  of  combat  days. 

COSAGE  output  data  has  been  made  available  by  the  Center  for  Army  Analysis 
(CAA)  as  (i)  the  total  forces  of  each  side,  Red  and  Blue,  by  platform  or  vehicle  type, 
present  at  time  0  for  each  of  several  specified  postures  (in  our  specific  numerical 
examples  we  analyze  Blue  Attack  vs.  Red  Hasty  Defense),  then  followed  (ii)  by  the 


corresponding  forces  at  time  2h  =  2  days  =  48  hours  later.  Data  furnished  concerning 
engagement  rates  and  kills  by  various  platforms  and  weapons  on  each  side  allow  one  to 
compute  average  engagement/shooting  rates  and  estimated  kill  probabilities  for  the 
various  Red  platforms  vs.  those  of  Blue,  and  vice  versa,  given  the  nature  of  the  platform- 
weapon-target  relation:  the  distinction  between  direct  “aimed”  fire  and  general  “area” 
fire.  Targets  are  in  principle  vulnerable  to  both. 

Several  new  discrete-time  dynamic  analytical/mathematical  models  have  been  fitted 
to  the  above  data  by  determining/estimating  simple  shooting  and  survival  parameters 
from  the  above  data.  These  simple  parametric  representations,  so-called  meta-models,  are 
then  used  to  compute/estimate,  in  spread-sheet  format,  future  force  sizes  and 
compositions  that  result  from  mutual  attrition,  as  well  as  the  corresponding  expenditures 
of  ammunition.  No  explicit  C4/ISR  or  Information  Operations  modeling  has  been  carried 
out,  but  such  effects  as  are  present  in  COSAGE  simulations  are  implicitly  represented  in 
the  parametric  representations  proposed.  It  is  then  shown,  for  instance,  that  tactical-level 
shooting  doctrine  can  be  modified  with  trivial  changes  in  kill  probability  parameters  in 
the  meta-models:  changes  from  kill  probabilities  associated  with  single  shots  at  acquired 
targets,  to  salvoes,  to  shoot-look-shoot,  etc.,  can  be  accomplished  off-line  and  then 
numerically  inserted  into  the  meta-model  computations.  This  opportunity  renders  painless 
certain  kinds  of  tactical  what-if  analyses.  For  more  detailed  discussion  see  Appendix  A, 
Firing  Options  and  Immediate  ( Imperfect )  BDA 

The  meta-models  are  also  capable  of  economically  and  efficiently  representing  likely 
shooting-rate  re-allocation  to  account  for  target  attrition.  The  adjustments  made  are 
speculative,  or  else  related  to  that  in  ATCAL,  an  Army  model  analogous  to  those  we 
propose.  We  provide  some  numerical  illustrations  of  these. 

We  can  also  represent  dynamic  feedback  adjustment  of  fire  rates  that,  if  overdone, 
could  result  in  unstable,  even  chaotic,  behavior. 
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Fitting  analytical  meta-models  (DISC-O-TIC)  and  examination  of  the  results  as 
compared  to  the  COSAGE  data  for  the  initial  2-day  calibration  period  can  be  instructive. 
For  example,  in  the  current  posture  of  Blue  Attack  and  Red  Hasty  Defense  the  given 
meta-models  for  Blue  survivorship  all  tend  to  fit  the  COSAGE  data  better  than  do  those 
for  Red.  This  is  likely  to  be  the  result  of  far  greater  attrition  of  Red  platforms  than  Blue 
(relatively  few  Red  survivors);  small  numbers  of  remaining  targets  tend  to  induce  extra 
shooting  unless  preventive  steps  are  taken;  this  effect  is  noticeable  on  a  suitable  graph.  A 
second  example  has  been  noticed,  in  which  COSAGE  attrition  was  significantly  greater 
than  that  of  the  fitted  model,  is  that  of  a  Blue  tracked,  lightly-armored  forward  scouting 
vehicle,  the  UM3CFV.  The  considerable  vulnerability  of  this  platform  stands  out  when  a 
platform-by-platform  plot  of  COSAGE  vs.  meta-modeled  survivorship  is  examined. 

It  is  strongly  recommended  that  some  form  of  meta-modeling  be  used  to  back  up 
every  complicated  simulation.  Taking  the  results  of  complicated  simulations  at  face  value 
without  intensive  diagnosis  and  critical  scrutiny  can  be  seriously  misleading,  particularly 
if  later  users  rely  on  the  report  that  the  tool  has  previously  been  used  by  others  and 
“therefore”  need  not  be  critically  examined  for  a  new  application. 
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1.  Introduction  and  Motivation 

Under  various  circumstances  it  is  natural  to  model  combat  attrition  in  discrete  time, 
i.e.  in  equal  time  steps  of  1-day  (24  hours)  or  even  of  2  or  3-day  durations.  This  is 
because  available  data  from  higher-resolution  tools,  such  as  the  Army’s  COSAGE,  may 
be  reported  on  a  daily  or  less  frequent  basis,  and  also  because  realistically-sized  model 
studies,  e.g.  of  theater-level  campaign  scale  may  require  iterative  computations 
(optimizations,  or  goal  programming)  that  are  best  performed  on  a  fixed  cycle  schedule, 
where  the  time  steps  are  not  too  small.  Attrition  (or  suppression)  calculations  must  fit 
smoothly  into  the  overall  model  computation,  which  is  often  on  a  coarse  time  scale. 

For  example,  work  being  earned  out  for  J-8  to  quantify  ammunition  consumption  in 
various  postures  involves  a  variety  of  computational  operations,  including  different 
attrition  calculations,  as  time  advances.  Consequently  it  is  efficient  to  allow  initially  for 
rather  sizable  time  steps,  on  the  order  of  several  days,  to  keep  computational  effort  under 
control. 

Traditional  combat/attrition  modeling  has  overwhelmingly  often  been  carried  out  in 
“continuous  time,”  often  using  the  language  of  differential  equations;  see  Taylor  (1983). 
The  best  known  of  these  are  the  classical  Lanchester  equations,  which  have  been  vastly 
generalized  (and  roundly  condemned,  with  some  justification).  Of  course  all  but  the 
simplest  combat  situations  prohibit  use  of  “closed  form”  solutions  to  such  equations,  so 
computational  solutions  have  been  developed;  these  actually  are  carried  out  by 
discretizing  time,  often  on  a  very  fine  scale.  Numerical  studies  that  involve  accurate 
solution  of  such  equations  can  be  quite  time  consuming,  especially  if  they  are  only  part  of 
a  much  bigger  effort. 

Both  data  availability  and  computational  burden  thus  urge  consideration  of  discrete¬ 
time  (daily,  or  several-day  time  step)  model  formulations  for  the  present  purpose.  An 
additional  argument  for  such  as  these  is  that  it  is  easier  to  enhance  the  simplest  (discrete 
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time)  models  by  introducing  adaptive  feedback  from  system  states  in  the  discrete-time 
format,  and  also  easier  to  accommodate  certain  stochastic  effects,  such  as  inter-agent 
visibility,  environmental  and  terrain  influences,  etc.,  than  it  is  in  continuous  time  (which 
models  must  be  discretized  to  compute  with,  in  any  case,  as  pointed  out  above).  We  carry 
out  some  such  enhancements  in  what  follows,  leaving  others  for  future  attention. 

Finally,  discrete-time  models  are  often  more  easily  understood  by  an  audience 
uncomfortable  with  the  ideas  of  calculus:  derivatives  and  integrals. 

2.  Model  1.1 

We  work  in  terms  of  state  variables  that  count/enumerate  the  numbers  of  Red  and 
Blue  opponents  of  various  types  and  capabilities  in  each  of  several  (sub)regions  at  times 
0,  h,  2h,  3h, 37  h, ....  However,  initially  we  discuss  the  simplest  such  model,  one  in 
which  a  single-type  Blue  force  is  in  combat  with  a  single-type  Red  force.  A  simplified 
deterministic  “fluid  approximation”  presentation  is  given  here,  and  made  general  in  later 
sections.  A  stochastic  version  is  presented  in  Appendix  A. 

2.1  Stripped-Down  Illustration:  Discrete-Time  Direct/Aimed/Allocated  Fire 

Begin  by  assuming  that  Red  (enemy)  and  Blue  (own)  forces  confront  each  other  in  a 
region  with  unspecified  dimensions  (but  only  temporarily;  generalization  follows).  The 
forces  are  of  the  same  type  here,  e.g.  are  both  tank  forces.  During  equal  consecutive  non¬ 
overlapping  time  periods  of  duration  h  (e.g.  single  day,  or  consecutive  several-day 
periods)  elements  or  the  forces  search  for  each  other  and  exchange  shots.  Assume  an 
aimed/assigned  fire  discipline  on  both  sides:  at  time  t  (e.g.  day,  3-day  step)  the  Blue  force 
size  B(t)  is  effectively  divided  by  the  number  of  Reds,  R(t)  (both  {B(t)}  and  {R(t)}  are 
treated  here  as  deterministic  fluids  -  the  “mean  field  approximation”  that  can  reasonably 
approximate  the  mean  of  a  Markov  chain  when  numbers  are  large).  Let  B(t)/R(t)  be  the 
number  of  Blues  in  a  specified  group  that  are  assigned  to  each  particular  Red  at  the 


5 


beginning  of  period  t  (fractional  numbers  are  admitted  and  interpretable  by  assuming  a 
shooter  processor-sharing  strategy).  Thus  each  Blue  has  his  exclusive  target  list  in  period 
t\  re-assignments  of  live  Shooters  (Blue  or  Red)  take  place  each  period.  On  average,  each 
Red  gets  the  attention  of  -  B(t)/R(t)  Blues  during  time  period  (t,  t  +  h );  each  Red  thus  may 
experience  B(t)/R(t)  •  nBR(*)h  shots,  where  nBR{9)h  denotes  the  rate  of  shot-fire  by  each 
member  of  a  Blue  subgroup,  the  notation  refers  to  the  possibility  of  describing  the 
effects  of  other  influential  variables,  such  as  the  number  of  live  Red  targets:  if  that 
number  becomes  small  then  the  rate  of  fire  at  it  might  also  be  reduced.(as  in  ATCAL:  see 
Appendix  D).  If  nBR{*)  is  taken  to  be  the  (constant)  maximum  acquisition  and  shooting 
rate  of  a  Blue,  then  the  effective  Blue  shooting  rate  may  potentially  be  unrealistically 
large,  meaning  that  many  Shooters  (Blue)  simultaneously  “pile  on”  a  possibly  hapless 
victim  (if  the  “victim”  is  already  dead,  or  is  a  false  target  or  decoy,  the  Shooters  are 
hapless,  at  least  as  represented  by  the  first  model  of  this  type).  Of  course  this  is  a  possible 
tactical  option  that  is  worth  evaluation.  It  may  be  more  munition-economical  than  the 
“pile-on”  option. 

A  complementaiy  setup  is  assumed  for  Reds:  each  Blue  is  assigned  a  specific  group 
of  R(t)/B(t )  Red  targets  at  the  beginning  of  period  t;  the  specific  groups  are  disjoint,  as 
before,  with  R(t)/B(t)  assigned  to  each  Blue.  This  symmetry  need  not  be  the  case,  but  is 
assumed  here:  available  data  from  simulations  seen  by  these  authors  can  support  no  more 
detail.  What-if  exploration  using  the  meta-model  is  relatively  easy. 

Suppose  combat  occurs:  compute  as  though  each  side  fires  simultaneously  (a 
convenient  assumption,  given  the  basic  data  currently  available).  Postulate  that  a  single 
representative  kill  probability  prevails  for  Red  shooter  vs.  Blue  targets,  denoted  KrB,  and 
Krr  for  Blue  shooter  against  Red  target.  Let  the  corresponding  one-shot  survival 
probabilities  be  denoted  by  krb  =  1  -  KrB  and  kbr  =  1  -  kBr.  These  initial  assumptions 
will  be  relaxed  subsequently  to  account  for  other  sources  of  variability,  in  general 
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associated  with  the  changing  environment  and  both  systematic  and  random/stochastic 
many-source  effects.  However,  it  is  currently  impossible  to  assign  “hard”  numbers  to 
these  effects  since  data  are  not  available. 

The  probability  that  a  single  Red  survives  for  one  period  is  0/^(0 .  corre¬ 

spondingly,  the  probability  that  a  single  Blue  survives  is  K^(*)ft/?(0/s(') 

Thus  the  expected  number  of  Reds  that  survive  a  time  period  (t,  t  +  h)  is  the  number 
at  t  +  h,  so  these  equations  result 

(2.1) 
(2.2) 

(again,  no  reinforcements  or  withdrawals  are  modeled).  Given  the  parameters,  equations 
(2.1)  and  (2.2)  are  easily  solved  recursively,  e.g.  on  a  spreadsheet,  given  the  parameters 
Kbr  and  krb,  nBR(*)  and  nRB(*)  and  initial  forces  5(0)  and  5(0). 

The  above  roughly  represents  dynamic  mutual  attrition  between  two  homogeneous 
sectorized  or  “aimed-fire”  forces  in  a  fixed  “basic”  time  period,  e.g.  one  day  under 
coordinated  aimed-fire  conditions.  It  does  not  attempt  to  model  the  moment-to-moment 
progress  of  a  battle,  but,  given  daily  (say)  data  5(0),  R(h),  R(2h),  R(3h), ...;  5(0),  B(h), ... 
5(13/?), ...  one  can  deduce  numerical  values  of  the  (daily  or  other  period)  average 
parameters  kbr,  krb,  nBR(»)  and  nRB{»),  where  the  latter  are  treated  as  constants;  other 
options  are  candidates  for  exploration,  but  supporting  data  is  non-existent.  However,  with 
enhanced  models,  plausible,  if  speculative,  doctrine  may  be  investigated,  and  sensitivity 
testing  carried  out.  The  “data”  {R(jh),  B(jh)}  can  come  from  runs  of  much  higher 
resolution  models  (e.g.  COSAGE);  some  results  of  using  such  to  “fit”  the  models  are 
reported  subsequently.  The  “fitted”  models  are  subsequently  used  to  extrapolate  the  high- 
resolution  results. 


R{t+h)=R{t)rBf')hB{,vR{,) 

B{t+h)  =  B(t)JcnRf')hRWB^ 
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See  Appendix  A  for  a _(one  of  many  possible )  stochastic  version  of  the  above  model. 


2.2  Connection  with  Classical  Lanchesterian  Differential  Equations 

Suppose  that  the  time  step  duration,  h—$ 0.  Expand  the  exponential  in  two  Taylor 
terms,  letting  icBR  =  \-kbr 


R{t  +  h)  =  R(t) 


1  +  nBR  (•)  In  kbr  h  + 0(/z2) 


so,  subtract  R(t)  and  divide  by  h  to  find 

R{t+h)-R(t) 


=  nBR(*)lntcBRB(t)  + 


O (*!) 


and  in  the  limit  as  h  tends  to  zero, 


dR(t )  .  .  _  .  , 

-^  =  «*«(•)  In 


when  B(t).  is  positive;  otherwise  the  derivative  must  be  taken  to  be  zero,  an  important 
boundary  condition.. 

For  small  kbr,  lnKBR  «  -kbr,  and 


and  likewise 


dR(t) 

dt 


*-BRnBR 


mt) 


(2.3,a) 


dB(t ) 
dt 


RBnRB 


m^) 


(2.3, b) 


again,  when  the  rhs  is  positive.  The  equations  (2.3  a&b)  are  recognized  as  the  classical 
square-law  Lanchester  equations.  It  must  be  recalled  that  unless  the  shooting  rates 
"**(•)  and  «**(•)  are  explicitly  forced  to  depend  on  R(t)  in  (2.3,a)  on  B(t )  in  (2.3,b),  with 

the  right-hand  sides  set  equal  to  zero  when  /?(/)  =  0,  or  B(t)  —  0,  the  equations  are 
incomplete  and  do  not  respect  physically  appropriate  boundary  conditions;  without 
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imposing  these  constraints  “solutions”  can  go  negative  or  exceed  initial  force  levels!  Such 
conditions  are  automatically  respected  in  (2.1)  and  (2.2)  even  if  nBR(»)  and  are 

taken  as  (often  unrealistically)  constant  in  (2.1)  and  (2.2)  for  large  B(t)  divided  by  small 
R(t). 


2.3  A  Queuing  Interpretation  and  Generalized  Lanchesterian  Results 

It  may  be  reasonable  to  conceptualize  opponents  as  virtually  “queuing  up”  for 
detection  and  attack  by  each  other.  For  the  moment,  consider  a  continuous-time  ( h  —>  0) 
“fluid”  model  of  such  queuing.  For  instance,  think  of  R(t),  the  surviving  Reds,  as  being  in 
a  virtual  waiting  (and  service)  line  for  B(t)  Blue  servers.  Then,  if  one  uses  an  analytical 
approximation  to  queuing  delay  effect  introduced  by  Rider  (1976),  and  Agnew  (1976), 
and  later  discussed  by  Filipiak  (1988), 


maximum 
service  rate 
(Blue) 


«(>) 

t  +  R(l) 

queue  +  in-service 
(Red) 


(2.4, a) 


and  likewise 


dt  ' — v— l+£(r) 

maximum  v.  v  7 

service  rate  queue  +  in-service 
(Red)  (Blue) 

(2.4,b) 

Such  equations  also  arise  in  biomathematics,  e.g.  in  enzyme  kinetics,  (also  a  combat 

situation),  and  may  take  the  more  general  form 

<«(»).•  £ 

dt  l  +  R(t)/KR 

(2.5, a) 

dB(t)  c  R(t)B{t) 

dt  '*Rb  l  +  B(t)/KB 

(2.5, b) 
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where  KR  and  KB  are  Michaelis-Menten  constants;  see  Murray  (1989).  Their  presence 
allows  further  opportunity  for  fitting  to  experimental  (simulation)  data. 

The  virtue  of  this  type  of  formulation  is  that  it  gracefully  bridges  the  gap  between  a 
target-rich  environment,  e.g.  for  B  as  in  (2.4,a),  when  R  is  large : 


M>)  _ 

dt 


(2.6) 


which  yields  the  Lanchester  Square  Law  of  (2.3)  above;  if  R  is  small,  so  Blue  is  target- 
poor,  i.e.  Red  targets  are  much  less  easily  found  and  killed,  we  find  (roughly) 

,)  (2.7) 

which  is  the  Lanchester  Linear  Law;  here  the  appearance  of  the  product,  B(t)R(t), 
represents  the  random  contact  rate  between  opposing  force  elements.  Owing  to  the 
smooth  saturation  of  the  shooter,  the  search-kill-rate  £*Br  <  %br  in  the  former  equation  is 

smaller  because  more  time  is  spent  searching  for  and  identifying  targets  before  shooting. 
Of  course  the  same  behavior  may  occur  for  Blue;  but  there  may  well  be  asymmetry 
between  force  behaviors. 

Now  go  in  the  time-honored  Lanchester  direction.  Division  of,  say,  (2.4,a)  by  (2.4,b) 
gives 


dR  Ur  (1  +  B) 

dB  £rb  (1  +  R) 


which,  upon  integrating,  yields 


(2.9) 


a  naturally  blended  version  of  the  Square  and  Linear  laws.  The  Michaelis-Menten 
generalization  of  (2.5)  is  seen  to  produce 
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(2.10) 


which  allows  further  tuning  between  the  “pure”  laws.  Generalizations  are  possible,  see 
Gaver  and  Jacobs  (1997),  but  are  omitted  for  the  present.  Taylor’s  books  (1980,  1983)  are 
basic  sources  for  this  topic. 


2.4  Extra/Over-Variability  in  Model  Response 

The  previous  model  types  are  not  now,  but  can  be  made,  explicitly  responsive  to 
either  systematic  and  explanatory/regression  variables,  or  to  random  (“hidden”)  sources 
of  variability  within  a  scenario  or  posture.  Among  such  system  effectors  can  be  range  and 
visibility  conditions  between  Shooter  and  Target,  and  maneuvers  taken  to  change  these. 
One  approach  to  include  these  effects  is  to  introduce  a  (one  or  more)  parsimonious 
parametric  modifications  of  the  survivor  probability,  k  ;  the  latter  may  be  patterned  on 
the  medical -biological  survivor-probability  methodology,  Cox  and  Oakes  (1984),  that 
incorporates  a  hazard  function,  deterministic  and/or  random.  This  leads  to  introduction  of 
a  revised  survival  probability 

*=(?)  =  g-*M0*(0)  (2.H) 

here  h(t;  x,  £ (t))  is  a  (partially)  random  hazard  rate.  In  classical  military  language  this 
resembles  (but  is  complementary  to  and  generalizes)  a  Killer-Victim  Scoreboard-,  it  is  a 
generalization  that  includes  systematic  observable  explanatory  variables,  plus  over¬ 
variability  from  many  causes. 

Example  1:  Range/Distance  Effects  Explicitly  Modeled 

The  data  available  to  parameterize  the  type  of  model  discussed,  e.g.  COSAGE  output, 
may  well 
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(a)  record  dBR(t),  a  characteristic  range  between  (segments  of)  the  two  forces  (here 
considered)  at  time  t  =  kh.  Other  things  being  equal  (they  will  not  be!)  we  anticipate  that  a 
hazard  component  due  to  range  alone  can  have  the  form 

KBR{t)  =  KBRf{dBR)  (2.12) 

where  the  basic  constant  kbr  and  the  function,  /(•),  of  distance  are  parameters  to  be 

determined.  This  form  allows  the  observed  survival  probability  to  be  exactly  matched  at 
two  specific  ranges,  and  be  adjustably  decreasing  with  increasing  range,  as  dBR(t) 
increases.  Of  course  other  explanatory  variables  such  as  terrain  cover  and  its  usage  also 
influence  kill  probability.  Unfortunately,  the  information  available  from  COSAGE  at 
present  kill  prevent  explicit  representation  of  such  effects. 


Example  2:  Random  Period-to-Period  Variation 

Think  of  R(t,  £ (r))  as  the  expected  or  mean  value  of  Red  force  size,  conditional  on 
£(t),  a  random  environmental  effect  that  persists  throughout  the  /th  period;  for  the  moment 
take  (f(r),  /  =  h,  2h, ...,  }  to  be  independently  and  identically  distributed  (iid);  however, 
the  analytical/mathematical  difficulties  do  not  increase  at  all  if  the  distributions  of  the 
environmental  effects  are  not  identical,  e.g.  depend  on  time  but  remain  independent,  //"the 
£(t)s  are  common  to  B(t)  and  R(t)  throughout,  as  could  represent  common  weather 
conditions,  then  we  can  put,  conditionally, 

/?(M(1))  =  R(0)Kl(x)nBR{')hB{0)/m 

(2.13) 

5(/z,f(l))  =  B(0)Jc£R{z)nBR{')hm/m 

Remove  the  condition  on  f(l)  to  obtain 
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R(h)  =  R(0)Ee  ee(')n*«(')h(m/m)hKB*  j  s  (.)fc(fl(o) /R(0)) In  kbr ] 

r  ,  (2-14) 

B(h)  =  5(0)££[/(1)',-(*)^(0)/B(°Wln^  ]  =  5(0)T[n,fi(.)/J(/?(0)/JB(0))ln^B] 
where  '?(•)  is  the  Laplace  transform  of  £(?), 

'F(*)  =  Ee[<T“(,)] 

and  where,  in  the  present  application, 

•s  =  -nBR(*)h(B(t)/R(t))lnlcBR 
>0  since  lnK^cO. 

On  the  fluid-approximation  principle,  see  Appendix  A,  calculate 

R(,+h) --  = 

r  -  -  ,  (215) 

Note  that  the  above  only  focuses  on  the  mean  survivorship  of  each  side.  In  fact,  the 
commonality  of  the  £(r)-effect  induces  positive  association  or  “correlation”  between  the 
attrition  rates.Operational  implications  of  this  phenomenon  are  not  investigated  for  the 
present. 

Illustration:  Let  {£(?)}  be  iid  and  Gamma,  with  mean  1  and  shape  parameter  /?, 
0  <  fi  <  °o,  (variance  1//?).  Then  since 

<2-i6) 

we  obtain 

R{t  +  h)  =  R(t)  1-”^  (R(t)/B{t)lnlcBR)  (2.17) 
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(2.18) 


B(t  +  h)  =  B(t) 


lnrM) 


-P 


Numerical  results  are  easily  obtained  using  spreadsheet  techniques.  Statistical  analysis  of 

data  to  support  this  model,  i.e.  to  estimate  the  value  of  beta,  has  not  yet  been  undertaken. 

However,  the  model  can  be  used  as  a  what-if  tool  to  study  the  numerical  effect  of  a 

plausible  model  change  on  survivability  and  munitions  expenditure. 

The  dependence  of  the  (mean)  survivorship  of  Red  and  Blue  is  seen  to  be 

systematically  affected  by  the  dispersion  of  the  ^-distribution.  The  mean  of  that 

distribution  has  been  fixed  at  unity.  If  /?  becomes  large  (/?-+  oo)  the  variance  (1//?)  hence 
standard  deviation  (l /V/?)  of  £(0  becomes  negligible,  and  the  original  model  reappears. 


If  P  is  small,  i.e.  close  to  or  below  unity,  the  survivorship  shape  changes,  and  radically  if 
P<  1:  relatively  high  survivorship  occurs  for  small  effective  values  of  killing  rate 
nBR{9)\nKBR ,  with  relatively  slower  trail-off  as  -nBR{*)lnK BR  increases. 

The  above  model  can  be  generalized  in  various  ways.  For  example,  the  £(/)-effect 
applied  to  one  side  for  several  specified  periods  can  represent  the  use  of  an  obscurant  or 
gas  (necessitating  use  of  movement-inhibiting  protection)  during  several  days  of  the 
campaign.  This  certainly  tends  to  reduce  the  capability  of  the  Shooters  on  the  receiving 
side  of  the  obscurant.  The  effectiveness  can  be  quantified  in  the  present  simplified  model, 
but  is  worth  doing  only  if  a  realistic  variety  of  weapon  types  is  in  operation.  This  more 
realistic  situation  is  a  candidate  for  consideration  in  later  work. . 


3.  Model  2:  Multitype  Shooters  and  Targets 

We  work  in  terms  of  state  variables  that  count/enumerate  the  numbers  of  Red  and 
Blue  opponents  of  various  types  and  capabilities  at  times  0,  h,  2  h,  3  h, 37  h, ... 

Note:  high-resolution  data  by  location  in  a  (sub)region  are  often  not  reported  as  output  in 
existing  simulation  models  (such  as  COSAGE),  although  the  model  itself  evolves 
spatially.  Analytical  opportunities  are  lost  by  hiding  the  spatial-temporal  aspect  of 
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combat  from  the  COSAGE  user.  It  is  recommended  that  more  such  detail  be  made 
available  to  the  user/analyst. 

3.1  State  Variables 

RQr',  t )  =  Number  of  Red  platforms  of  type  jR  in  region  iR  at  time  t 
B(jB;  t )  =  Number  of  Blue  platforms  of  type  jB  in  region  iR  at  time  t 

3.2  Parameters 

Pb(/b,  vvb;  t)  =  rate  of  fire  of  weapon  type  wB  from  Blue  platform  type  jB  in 

region  iB  during  it,  t  +  h ] 

MbrUb,  wb\  jp,  t)  =  fraction  of  weapons  of  type  wB  from  platform  of  type  jR  fired  at 

Red  platforms  of  type  jR  during  the  time  interval  [t,  t  +  h]. 

Kbr(Jb,  wb\  Jr)  =  probability  of  kill  for  Blue  platform  of  type  jB  firing  with 

weapon  wB  at  Red  platform  of  type  jR 

Note:  Attrition  need  not  be  permanent  (“kill”):  by  expanding  the  state  space  we 
accommodate  partial  (e.g.  mobility)  kills,  and  perhaps  temporary  psychological  kills 
(“suppression”);  however,  recovery  from  suppression  is  not  explicitly  modeled  here,  but 
to  do  so  is  not  difficult.  Again,  data  support  for  this  is  not  available  explicitly  from 
current  COSAGE  data. 

3.3  Aimed  or  Direct  Fire 

This  model  type  is  useful  for  representing  the  attrition  of  Red  by  Blue  (and  vice  versa) 
when  Blue  weapons  can  be  allocated  to  exercise  “aimed”  or  “direct”  fire  at  batches  of 
Reds  whose  type  and  location  is  presumed  known  (visible  either  directly  or  indirectly)  to 
Blue  shooters.  Tank  warfare  is  an  example.  The  total  shots  capable  of  being  fired  at  Red 
targets  of  type  7^  by  Blue  platforms  of  type  jB  firing  weapon  wB  in  (t,  t  +  h]  is 

nbr  (  jB  ,wB-jR;t)  =  B(jB ;  t  )pB(jB ,  wB ;  t)d)BR(jB ,  wB ;  jR ;  t)h 

(3.1) 

=  B(jB  ;t)[nBR(jB ,  wB ;  jR ;  t)h] 
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In  words,  nBR(*)h  is  the  number  of  shots  capable  of  being  fired  by  an  individual  type  jB 
Blue  platform  at  t  against  Reds  of  type  jR  during  (t,  t  +  h\.  The  number  of  Reds  of  type  jR 
at  time  t  +  h  (no  reinforcements)  is  obtained  by  the  argument  of  Section  2  preceding:  the 
basic  notion  of  aimed  or  direct  fire  is,  here,  that  a  certain  number  of  the  total  shots  fired 
by  a  platform  in  a  given  category  is  allocated  to  each  single  Red  unit,  type  jR.  Only  a 
subset  of  Blues  is  allocated  to  each  of  a  particular  subset  of  Red  (Blues  do  not  target 
totally  at  random,  nor  do  Reds).  In  the  present  deterministic/fluid-like  approximation  this 
suggests  that  each  one  such  Red  is  the  target  of  a  potential  NB(js,  wB;  jR,  t)/R(jR,  t )  shots. 
The  probability  that  such  a  target  survives  for  one  /z-period  is 

JB  ™B 

(3.2) 

JB  WB 

from  which  come  the  survival  recurrences  (assuming  both  sides  behave  in  the  same  way). 

Model  2.1: 

(3. 3, a) 
(3.3, b) 

where  the  probability  of  a  Red  of  type  jR  surviving  aimed  fire  (AF)  for  one  time  unit,  h,  is 
thus  given  by  the  right-hand  side  of  (3.2).  An  equivalent  formula  holds  for  Red  DF  vs. 
Blue.  We  now  abbreviate  the  number  of  shots  per  Blue  at  a  Red  jR  as  nB(*,  jB,  wB;  jR;  t). 
As  before,  in  Section  2,  we  use  the  to  indicate  an  unspecified  (to-be-specified,  but 
otherwise  taken  to  be  constant,  usually  representative  of  the  first  few  days  of  COSAGE 
data)  dependence  of  the  number  of  shots  on  either  deterministic  or  stochastic  sources  of 
variation. 


B(jB  ;t+h)  =  B(jB ;  t  )S(jB ;  t) 
R{jR‘,t +h )  =  R(jR;t)S(jR;t) 
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3.4  Firing  Rate  Allocation 

As  remarked  in  Section  2,  it  is  not  satisfactory  or  realistic  to  let  a  firing  rate,  e.g.  of 
Blue  weapon  wB  from  a  platform  type  jB,  against  a  Red  of  type  jR,  i.e. 

nBR{jB,wB;jR,t ) 

remain  constant  over  time,  as  opponent  numbers  change  {decrease  only  in  this  model):  if 
B(jB,  t)  becomes  relatively  large  compared  to  R(jB;  t)  then  an  unrealistic  number  of  shots 
may  be  fired  at  a  Red  target  —  one  that  in  practice  is  either  already  dead,  or  may  have 
moved! 

The  numbers  of  shots  of  various  types  fired  at  a  perceived-alive  target  are  actually 
decision  variables.  The  decision  is  guided  by  rules  of  doctrine,  and  training,  but  must  rely 
in  real  life  on  a  tank  commander’s  individual  skill  and  perceptiveness.  Our  present 
models  do  not  directly  account  for  variability  of  operator  performance.  In  the  context  of 
the  current  model  type,  and  others  more  detailed  and  high-resolution,  we  will  suggest  and 
test  some  rules  that  behave  with  qualitatively  correct  properties.  For  the  first  model  we 
simply  take  firing  rates  as  constants,  estimated  from  data  (COSAGE).  This  assumption  is 
subsequently  modified. 

Some  Shot  Allocation  Rules 

The  “historical”  COSAGE  run  record,  of  2h=48  hours’  duration  for  this  report’s 
sample  data  set  has  been  analyzed  to  develop  estimates  of  shooting  rates  and  kill 
probabilities  for  platform- weapon  combinations  against  specified  target  types. 

As  noted,  in  the  first  model  the  above  rates  and  probabilities  are  not  altered  (they 
remain  constant)  over  time:  the  same  shooting  rates  per  specified  Blue  Shooter  vs.  Red 
target,  and  vice  versa  are  assumed  to  prevail  throughout:  these  are  the  Shooter  vs.  Target 
rates  that  characterize  the  initial  48  hours  of  the  COSAGE  simulation  data.  It  is  justifiable 
to  question  this  assumption,  since  presumably  total  rate  of  fire  directed  at  individual 
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targets  will  diminish  as  the  number  of  targets  per  assigned  shooters  decreases.  This  may 
occur  because  of  increased  difficulty  of  acquiring  alive  eligible  targets. 

Here  are  some  rules: 

(a)  Initial  (0,2h]  =  48  hour  shooting  rates  per  shooter  maintained  throughout  the  elapsed 
battle  time. 

(b)  Initial  shooting  rates  redistributed  each  time  period  across  alive  targets,  in  same 
proportion  as  during  initial  period  (0,  2h]. 

(c)  Initial  shooting  rates  redistributed  each  time  period  in  proportion  to  number  of 
surviving  targets  in  each  target  class.  The  redistribution  need  not  be  proportional  but 
can  be  some  general  or  weighted  function  that  is  selected  for  control  purposes. 

(d)  Initial  shooting  rates  redistributed  according  to  the  perceived  lethality  (against 
appropriate  associated  shooter  platforms)  of  each  target  class.  Approximately  and 
myopically:  weight  the  surviving  numbers  of  each  target  class  by  the  probability  of 
kill,  or  rate  of  kill,  against  its  shooter-class  targets.  In  short,  a  shooter  of  either  side 
may  well  tend  to  shoot  first  at  its  currently  most  dangerous  opponent. 

(e)  Initial  gross  shooting  rates  maintained  constant  in  each  subsequent  shooter-target 
pairing.  This  is  equivalent  to  applying  the  logic  of  the  ATCAL  model  assumptions  in 
each  time  period;  see  Chap.  6,  expression  (6.2.1)  of  Caldwell,  et  a\.  (2000).  In  turn, 
this  amounts  to  adjusting  the  engagement-killing  rate  so  as  to  keep  the  target  type 
attrition  rate  constant  during  each  day  of  the  battle,  which  leads  to  an  exponential 
decrease  in  targets  throughout  the  battle. 

Comments:  There  appears  to  be  no  basic  reason  why  such  a  rule  is  especially  natural 
or  “optimal”.  It  does  appear  to  be  part  of  one  standard  Army  analytical  tool,  namely 
ATCAL,  and  hence  is  worth  examination  in  the  present  context.  It  implies  a  certain 
shooter  restraint:  if  the  number  of  shooters  vs.  members  of  a  class  of  alive  targets 
becomes  greater  than  was  the  case  initially,  then  the  shooters’  shooting  rate  is  reduced  so 
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that  the  original  shooting-killing  rate  prevails  in  the  new  time  interval.  This  may  be 
qualitatively  plausible  in  some  rough  sense.  Why  the  adjustment  to  achieve  exponential 
decrease  of  target  forces  might  be  descriptive  of  real  combat  remains  mysterious. 

The  basic  direct-fire  allocation  of  Blue  platforms  to  Red  targets  in  the  present  model 
begins  by  assigning  a  set  of  Blue  platforms  exclusively  to  a  particular  set  of  Reds. 
Consequently,  not  every  Blue  platform  in  the  situation  represented  can  acquire  or  shoot  at 
any  Red  platform;  there  may  be  different  assignments  throughout  a  basic  period  of 
duration  h  (e.g.  a  day),  but  it  is  presently  assumed  that  those  assignments  are  respected. 
Given  the  data  available,  it  is  the  basic  current  assumption  that,  during  period  ( t ,  t+h) 
each  Red  of  type  jR  is  a  candidate  to  experience  shots  from  some  Blues  of  type  jB.  e.g. 
during  period  t  (e.g.  day  t  to  t+h)  each  Red  of  type  jR  receives  on  average  Blue  shots  of 
weapon  type  wB  from  platforms  of  type  jB  that  number 


„  .  ,,  .  y4n(»,jB,wB;jR;t)B(jB;t) 

Blue  wB  shots  at  a  Red  of  type  Jr  \  / 


(3.4) 


during  time  period  (t,t  +  h)  R(Jr  ;  t) 

Data  available  from  COSAGE  runs  provides  the  possibility  of  calculating/estimating  the 
initial  shooting  rate  of  a  Blue  type  jB  against  a  Red  target,  type  jR: 

#  of  wB  -  shots  by  all  B(  jB  ,0)  in  (0,h) 


^(•,jB,wB,jR;h)  = 


B{jB,  0) 


(3.5) 


The  rate  of  shooting  is  very  likely  to  be  strongly  dominated  by  the  rate  of  target 
acquisition:  the  time  between  shots  is  mainly  the  time  between  successful  detection, 
acquisition,  and  classification  (or  misclassification)  of  a  target;  the  actual  shooting  time 
(weapon  firing  and  delivery  on  target)  is  negligible  by  comparison. 


3.5  Poisson  Shooting-Rate  Model:  Model  2.2 

Invocation  of  the  Poisson  process,  Feller  (1968),  provides  a  plausible  and  useful  next- 
stage  model  that  incorporates  stochastic/chance  variability  into  target  acquisition  and 
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shooting.  During  a  generic  period  ( t ,  t  +  h)  a  group  of  Shooter  platforms  allocated  to  fire 
specified  weapons  at  a  specified  target  type  search  for  their  quarry  in  such  a  way  that  the 
number  they  find  is  random  (realistically),  and  is  governed  by  a  Poisson  random  process. 
If  the  Shooter  is  of  type  js  and  the  quarry  of  type  jg  (here  if  Q  =  B,  S  =  R,  while  otherwise 
Q  =  R,  S  =  B )  the  random  number  found,  per  Shooter  of  type  js  in  (t,  t  +  h )  has  the 
Poisson  distribution  with  rate  X sQ{*,js,jQ );  if  the  Poisson  random  version  of  nsQ(*,js,jo)h 
has  mean  XSQ(»,js,jQ)h  then  it  can  be  shown  that  (3.3)  is  replaced  by  (3.6,b)  below,  with 

-XX-M*  Jr Jb )h  KRB(jR ,wR  ;jB )(R(jR ,t)/B(js ,<)) 

bUb ,t  +  h)  =  B(jB , t)e  j*  w*  (3.6,a) 

and 

”X  **R  (*Jb  'wb  Jr  )*  *BR  (jB  .H-fi \jR )( B(jB  ,1  )/R(  jR ,/)) 

R{jR,t  +  h)  =  R(jR,t)e  jBwB  .  (3.6, b) 

Here  we  have  incorporated  the  fraction  of  different  weapons  shot  by  the  same  platforms 
directly  into  the  parameterized  search  rate.  As  before,  the  symbol  denotes  the  possible 
dependence  of  the  acquisition  latency  on  other  explanatory  variables,  such  as  location, 
range,  previous  acquisitions  experienced,  etc.  These  variables  are  not  included  here 
because  they  cannot  be  obtained  from  current  COSAGE  data. 

4.  Connection  with  COSAGE  Data 

Note  that  the  raw  output  of  COSAGE  data  currently  made  available  aggregates  Red 
(and  Blue)  force  types  over  their  various  locations. 

COSAGE  data  from  actual  runs,  or  other  historical  data,  may  give  numerical  values 
for  the  fraction,  cobr  of  weapons  of  type  wB  historically  fired  by  Blues,  type  jB,  vs.  Reds, 

type  and  vice  versa.  The  data  allows  empirical  determination  of  the  systems’  attrition 
power. 

An  alternative  is  to  (eventually)  simulate  rule-based  adaptivity:  the  rate  of  fire  at  a 
given  batch  of  Red  targets,  e.g.  tanks  at  a  certain  location,  may  be  made  much  higher  than 
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at  APCs  at  another  location;  this  can  be  achieved  by  adjusting  the  functional  form  of 
C0BR(jB,  wb-Jr)-,  the  latter  can,  and  should,  depend  on  (currently  unmodeled)  availability  of 
ammunition  in  various  categories,  and  of  course  the  presumed  threat  of  the  different  Red 
target  batches. 

4.1  Adapting  the  Model  to  Easily  Utilize  COSAGE  Data 

Models  used  to  estimate  the  expected  number  killed  in  direct  fire  using  parameter 
estimates  from  COSAGE  are  summarized  below. 

Notation 

R(Jr,  t)  =  Number  of  Red  shooters  of  type  jR  at  time  t 
B(j'b,  t)  =  Number  of  Blue  targets  of  type  js  at  time  t 

NRB(jR,  wr-Jb,  t)  =  (Mean)  number  of  shots  fired  by  all  Red  shooters  of  type  jR  using 
munition  wR  at  all  Blue  targets  of  type  jB  during  period  (0,1] 

nbr(j'b,  wb-Jr,  t)  =  (Mean)  number  of  shots  fired  by  all  Blue  shooters  of  type  jB  using 
munition  wr  at  all  Red  targets  of  type  Jr  during  period  t:  (0,1] 

2 h  =  number  of  days  of  combat  represented  in  the  COSAGE  run 

r  =  the  number  of  COSAGE  replications;  often  16. 

nRB(jR,  wR;  jB,  t)  =  NrbUr,  wr;  jB,  2 h)/R(jR;  0)2 h  =ARB(jR,wR;  jB ) 

"brUb,  WB; Jr,  t)  =  NbrUb,  WB;  jR,  2h)/B(jB;  0)2 h  =ABR(jB,wB;  jR ) 

KrbUr,  wr-Jb,  t)  =  (Mean)  number  of  Blue  targets  of  type  jB  killed  by  shots  fired  by 
all  Red  shooters  of  type  jR  using  munition  w/?  during  (0,1] 

KbrUb,  wB;jR,  t)  =  (Mean)  number  of  Red  targets  of  type  jR  killed  by  shots  fired  by 
all  Blue  shooters  of  type  jB  using  munition  wB  during  (0,1] 

KRB(jR,wR;  jB,2h)+-  Bayesian  estimated  probability  of 

KRB(jR,wR;jB)  = - —L.=  kill  of  Blue  type  jB  by  Red  (4.1) 

N, rB  (jR ,  wR ;  jB  ,2 h)  +  -  platform,  type  jR  using  weapon 

r  type  w/j.  Uses  uniform  (0,1)  prior 
for  the  probability  of  kill 
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Kbr ( jB ,  wB ;  Jr  ,2 h) + -  Bayesian  estimated  probability  of 

kBrUb  ,wB;jR)  = - r-  =  kill  of  Red  type  jR  by  Blue  (4.2) 

Nbr  ( jB ,  wB ;  jR  ,2 h) + -  platform,  typ tjB  using  weapon 

r  type  wB  .  Uses  uniform  (0,1)  prior 
for  the  probability  of  kill 


First  model:  Model  2.1 

Rates  of  fire  observed  during  (0,2 h]  in  COSAGE  are  used  for  all  future  times.  The 
same  for  kill  probabilities.  Both  are  estimated  as  described  under  Notation  above. 

)R 


Poisson  model:  Model  2.2 

B{  jB ,  t  +  h)  =  B(  jB ,  exp{-/UB  ( jR ,  wR ;  jB  )h  krb  (  jR ,  wR ;  jB  )R(jR  ,t)/ B(jB,  t)}  (4.4) 

jR 

Expected  number  kills  model:  Model  2.3 

B(js  ,t  +  h)  =  B(jB,t)~Y,  nRB  ( jR ,  wR ;  jB  )R(jR ,  t)hKRB  (jR ,  wR ;  jB )  (4.5) 

Jr 

if  B(jB,t  +  h)> 0  and  0  otherwise. 

Fluid  queuing  model:  Model  2.4 


B(  j B ,  t  +  h)  =  B{  jB ,  t )  -  £  nRB  ( / R ,  wR ;  jB ) /?( j R ,  t) - — — 1 — r  hK RB  ( jR ,  wR ;  jB ) 

jR  1  +B{jB,t) 


if  B(jB ,  t  +  h)  >  0  and  0  otherwise. 

Random  Environment  Gamma  model:  Model  2.5 


(4.6) 


B(jB,t  +  h)  =  5(yB,r)|^[ 

jR 


1- 


nRB(jR ,  wR ;  jB )  ln(l  -  k  RB(jR  ,wR;jB  ))R(jR  ,t)h/  B(jB,t) 

~P 


Figures/graphs  to  be  presented  graphically  present  the  COSAGE  data  plotted  vs. 
the  fitted  models.  These  are  results  comparing  the  average  number  of  LIVE  platforms 
surviving  from  COSAGE  and  projected  by  the  above  models  for  each  platform.  Such 


(4.7) 
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graphs  are  useful  diagnostic  tools  that  can  reveal  anomalies  in  data  (and  point  out 
programming  errors!).  COSAGE  summary  output  reports  summary  statistics  averaged 
over  a  specified  number  of  replications.  For  aimed  fire  weapons  the  statistics  include  the 
average  number  of  weapons  of  type  w  fired  by  platform  of  type  j  at  a  target  of  type  k  and 
the  average  number  of  targets  of  type  k  killed  by  weapons  of  type  w  fired  by  platforms  of 
type  j  for  a  2-day  period.  Also  specified  are  the  initial  numbers  of  Red  and  Blue 
platforms,  the  number  of  replications  and  the  combat  posture.  The  kill  probability  of  a 
single  aimed  fire  weapon  wB  fired  by  a  Blue  platform  of  type  jB  against  a  jR  Red  in 

situation  j  can  be  estimated  from  the  COSAGE  output  as  follows. 

Let  KBR(jB,wB.jR;s)  be  the  average  (e.g.  over  16  replications)  number  of  Reds  of  type 

jit  killed  by  aimed  fire  weapon  wB  fired  by  Blue  platforms  of  type  jB  for  the  two-day 
period  during  posture  j.  Let  NBR(jB,wB;jR;s)  be  the  average  number  of  aimed  fire 

weapons  wB  fired  by  Blue  platforms  of  type  jB  at  Reds  of  type  jR  for  the  two-day  period 
during  posture  s.  In  aimed  fire  each  weapon  is  assumed  to  be  able  to  affect  exactly  1  Red 
target;  we  ignore  weapons  fired  at  already  killed  targets.  An  estimate  of  the  kill 
probability  of  a  single  aimed  fire  weapon,  type  wB,  shot  by  a  Blue  of  type  jB  against  a  jR 
Red  in  posture  5  during  the  two-day  COSAGE  run  is 


^BR{jB,wB;  jR\s) 


KBR(jB,wB-JR;s) 

NBR(jB,wB;jR;s) 


(4.8) 


(Note:  this  is  actually  the  maximum  likelihood  estimate  of  a  binomial  model’s  kill 
probability,  here  KBR(jB,  wB;jR;  s).  It  is  scenario-specific,  and  range  or  other  condition 
dependence  is  not  specified.) 

Now  to  avoid  having  to  estimate  a  kill  probability  as  0  owing  to  small-sample  bad 
luck  we  replace  Kbk/ Nbr  by  a  Bayes  estimate  that  assumes  kills  are  binomial  and  applies 
a  uniform  prior  (one  can  use  a  more  informative  prior  if  one  is  available);  the  result  is  the 
estimate  of  kill  for  Red  of  type  jR  when  fired  upon  directly  by  Blue  of  type  jB  using 
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weapon  type  wB  during  phase  5;  if  r  is  the  number  of  COSAGE  replications  (e.g.  16);  the 
Bayes  probability  of  kill  for  a  single  shot  is 


KBR{jB>WB>jR>s)  +  — 

KbrOb ,  ;  jR ; 5)  =  — - r- . 

A rBR{jB,wB;jR;s)  +  - 

r 

Thus,  the  probability  a  particular  Red  survives  a  shot  is 


A 

KbrUbWJr',  5)=  1- 


KBr  (jB,wB;jR;s)  +  — 
_ _ r_ 

_  2 

NBR(jB,wB;jR',s)-\ — 


(4.9) 


(4.10) 


4.2  Area  or  Unaimed  Fire:  Model  3.1 

Suppose  certain  weapons  fire  in  an  unaimed  area-fire  manner.  Then  (4.3)  must  be 
modified  to  reflect  the  way  that  any  shot  places  any  Red  target,  type  jR,  in  location  iR  at 
risk.  One  formula  will  handle  this  for  each  phase/stage  s;  put 

(4.11) 

JB  VVfi 

where 


[  1  if  weapon  B  is  aimed; 

0{wb)  =  \  (4.12) 

0  if  weapon  B  is  area. 

and  Nbr  is  defined  in  (3.1).  Of  course  the  corresponding  kill  probabilities  must  be  made 
weapon-specific.  It  is  even  possible  to  let  the  indicator  £(•)  depend  on  other  features  of 
combat,  such  as  distance,  general  terrain  features,  visibility,  etc.,  and  to  let  6{-)  take  on 
values  other  than  0  or  1.  There  is  currently  no  information  quoted  from  COSAGE  runs  to 
permit  such  refinements.  This  represents  a  lost  opportunity  for  learning. 
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4.3  Using  Available  COSAGE  Data  to  Estimate  Survival  for  Area  Fire  Weapons 

In  this  section  we  provide  formulas  for  the  effect  of  indirect  or  area  fire  analogous  to 
(3. 3, a),  (3.3,b),  (3.6,a),  and  (3.6, b)  for  aimed  fire.  We  use  only  data  now  available  from 
COSAGE  runs,  so  the  footprint-that-affects-several-targets  effect  must  be  inferred 
indirectly.  No  doubt  the  process  used  here  can  be  improved,  but  the  method  necessarily 
uses  data  presently  available  from  a  detailed  “realistic”  model,  COSAGE,  to  “fit”  a  much 
cruder  model:  simple  recurrences  (3.3, a)  and  (3.3, b). 

The  argument  below  applies  to  a  generic  area-fire  situation;  it  is  easily  made  specific. 

There  are  R  Red  targets,  all  of  the  same  type  (e.g.  tanks),  and  all  equally  vulnerable  to 
Blue  weapons.  Let  there  be  Sb  Blue  indirect-fire  weapon  shots  of  specified  type  aimed 
into  the  region  where  the  Reds  reside.  Let  or  be  the  probability  that  a  Red  (target)  is 
affected  (e.g.  susceptible  or  exposed  to  being  killed,  but  possibly  damaged  or  suppressed 
within  a  footprint)  by  a  Blue  indirect  weapon  shot.  Let  the  indicator  function  /,(/)  =  1  if 
the  1th  shot  affects  the  target  number  j;  otherwise  ///)  =  0.  Then  the  total  number  of  targets 
affected  by  indirect  weapons  during  a  replication  of  COSAGE  is  the  sum 

Hfitm.  <4 13> 

/=i  j= i 

a  random  variable,  the  observed  value  of  which  on  replication  k  is  n/k),  k  =  1,...,  r.  Now 
assuming  symmetry  over  the  area  fired  upon,  Ft///)]  =  a  for  each  shot  and  target 
(possibly  more  plausible  if  all  shots  are  fired  simultaneously  so  that  targets  do  not  move). 
Thus 


E[A/]  =  SsRor  (4.14) 

and  so  we  can  estimate  aby  moments: 

»/=-]£n/(*);  (4.15) 

r  k 
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consequently  an  estimate  of  a  is  a  =  nf/SBR;  the  number  of  shots  SB  is  estimated  by  the 

average  of  the  number  of  shots  observed  during  the  r  replications  of  48  hours  of 
COSAGE  combat. 

Suppose  £  is  the  probability  that  an  affected  (within  footprint)  target  is  killed.  The 
probability  that  a  particular  target  is  killed  is  thus  at,  (first  affected,  then  killed);  the 
probability  of  surviving  one  shot  is  the  1  -  at,  for  any  target  (symmetry).  Assume  now  that 
the  probability  that  a  single  target  survives  (all)  S  shots  is  (1  -  a£)s\  this  is  perhaps  more 
plausible  if  the  target  cannot,  or  does  not,  move.  In  that  respect  it  is  a  lower  bound  on  a 
particular  target’s  survival.  Put  Kw(a)  =  at,. 

COSAGE  records  the  average  number  of  targets  killed  by  the  indirect  fire  weapons 
over  the  r  (often  r—  16)  replications:  k/?(2)  for  Red,  while  R( 0)  is  the  initial  number  of 

Reds.  By  symmetry  and  the  method  of  moments  put 


kR(2)  =  R(0)SBKw(a)  (4.16) 

from  which  the  probability  that  any  Red  is  killed  by  a  single  Blue  shot  is  estimated  as 


v  '  r(o)sb 


(4.17) 


where  Sb  is  the  average  number  of  indirect-fire  (Blue)  weapons  shot  during  the  r 
replications. 

To  avoid  having  to  estimate  a  survival  probability  as  1  owing  to  small-sample  luck  we 
replace  k^(2)/ SBR( 0)  by  a  Bayes  estimate  that  assumes  kills  are  binomial  and  applies  a 
uniform  prior;  the  result  is  the  estimate  of  survival  for  Red  when  fired  upon  indirectly  by 
Blue,  using  weapon  type  wB:  if  r  is  the  number  of  replications  (e.g.  16), 


k  b/?(wb,  y‘/?)  = 


( 


1- 


fc«(2)  +  l/r 

R(0)5b  +2/7 


(4.18) 
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where  SB  =  B(0)nBR(-)  since  the  above  applies  to  Blue  shooting  at  Red.  A  symmetrical 
formula  holds  when  Red  is  indirectly  shooting  at  Blue. 

5.  Numerical  Illustrations 

We  provide  a  number  of  graphical  displays  of  COSAGE-data  supported  survivorship 
and  Munitions  expenditures  by  various  platforms.  Fitted  meta-models  are  compared  to 
COSAGE  data  for  the  two-day  time  period  for  which  COSAGE  data  is  available  for  the 
posture  Blue  attack-Red  hasty  defense.  Then  comparisons  between  model  types  are 
displayed.  The  model  types  considered  are  the  Classical  Survival  Model  1  (4.3),  the 
Poisson  Acquisition  model  (4.4)  and  an  ATCAL-like  attrition  model  (D.4)  and  (D.8). 
Comparison  of  Figures  A.2  and  B.2  show  that  the  Classical  Survival  Model  1  summarizes 
COSAGE  number  of  platforms  surviving  results  for  the  Blue  forces  better  than  for  the 
Red  Forces.  In  this  scenario,  the  Red  forces  suffer  greater  attrition  than  the  Blue  forces. 
Classical  Survival  Model  1  tends  to  underestimate  the  number  of  surviving  platforms. 
The  Blue  platform  with  the  largest  negative  discrepancy  in  Figure  A.2  is  Platform  25 
which  is  the  UM3CFV,  a  lightly  armored  tracked  forward  Blue  scouting  vehicle  that 
suffers  greater  attrition  in  the  Classical  Survival  Model  1  than  in  COSAGE.  Figures  C.l 
and  C.2  display  the  amount  of  Blue  munitions  used  and  the  number  of  Red  surviving 
platforms  after  8  days  for  the  Poisson  acquisition  model  versus  the  classical  survival 
model.  The  Poisson  acquisition  model  tends  to  predict  more  munitions  used  for  fewer 
Reds  killed.  Figures  D.  1  and  D.2  compare  the  amount  of  Blue  munitions  used  during  8 
days  for  classical  survival  model  1  and  the  Poisson  acquisition  model  with  and  without 
proportional  reallocation  of  fires  as  specified  in  equation  (E.l)  of  Appendix  E.  Not 
surprisingly,  the  re-allocation  of  fire  results  in  more  Blue  munitions  being  used.  However, 
the  effect  is  greater  for  Classical  survival  model  1  than  for  the  Poisson  acquisition  model. 
Figures  E.l  and  E.2  compare  the  number  of  surviving  Red  platforms  after  8  days  for  the 
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classical  survival  model  1  and  the  Poisson  acquisition  model.  As  expected  the  weapon 
reallocation  results  in  a  few  more  Red  platforms  being  attrited.  Once  again,  the  effect  is 
more  extreme  for  the  classical  survival  model  1.  Figures  F.l  and  F.2  display  the  number 
of  Blue  munitions  fired  during  8  days  and  the  number  of  Red  platforms  surviving  after  8 
days  for  the  classical  survival  model  1  and  the  ATCAL-like  model  described  in  Appendix 
D.  The  ATCAL-like  model  tends  to  predict  fewer  Blue  munitions  expended  during  the  8 
days  than  the  Classical  Survival  model  1.  The  ATCAL-like  model  tends  to  predict  fewer 
Red  platforms  surviving  for  those  Red  platforms  that  have  a  relatively  large  number 
surviving  in  Classical  Survival  model  1.  If  Classical  Survival  Model  1  predicts  a  small 
number  of  Red  platforms  surviving,  then  the  ATCAL-like  model  tends  to  predict  a 
somewhat  greater  number  of  platforms  surviving. 

As  they  are  currently  fit  to  the  particular  COSAGE-model  posture,  none  of  the  meta¬ 
models  differ  enormously  in  their  implications.  This  can  change,  for  example,  if  the  meta¬ 
models  are  numerically  parameterized  differently,  i.e.  if  the  Random  Environment  model 
(4.7)  is  used  with  the  value  P  small,  or  if  the  ATCAL  tuning  parameter  f  is  altered.  The 
models  require  exploration  that  shows  which  have  the  most  conservative  but  realistic 
assessment  of  kill  rate  for  munitions  expended. 
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Number  of  Blue  surviving  platforms  from 
classical  survival,  model  1 


A.l  Classical  Survival  Model  1  vs.  COSAGE  Survival,  for  Blue  (No 
reallocation  of  fire) 

Number  of  Blue  Surviving  Platforms  After  2  days 
Classical  Survival,  Model  1  versus  COSAGE  after  Attrition  by  Direct  Fire  Weapons 
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Fraction 


A.2  Classical  Survival  Model  1  vs.  COSAGE  Survival,  Normalized 
Difference  for  Blue 
(No  reallocation) 

Comparison  of  Numbers  of  Surviving  Blue  Platforms  for  Classical  Survival  Model  1 

and  COSAGE 

(COSAGE-[Classical  Survival,  Model  1])/(lnitial  Number  of  Platforms) 


Platform  Index 
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A.3  Classical  Survival  Model  1  vs.  COSAGE  Munitions  Expenditure:  Blue 

Number  of  Blue  Munitions  Fired  During  2  days 


B.l  Classical  Survival  Model  1  vs.  COSAGE  Survival,  for  Red 

Number  of  Surviving  Red  Platforms  after  2  days 
Classical  Survival,  Model  1  versus  COSAGE  after  Attrition  by  Direct  Fire  Weapons 
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Fraction 


B.2  COSAGE  Survival— Classical  Survival  Model  1:  Normalized  Difference 


for  Red 

Comparison  of  Numbers  of  Surviving  Red  Platforms  Classical  Survival,  Model  1  and 

COSAGE 

(COSAGE-[Classical  Survival  Model  1])/(lnitial  number  of  Platforms) 
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Number  of  munitions  used  Poisson  acquisition  model 


C.l  Munitions  Used  by  Blue  in  8  Days:  Classical  Survival  Model  1  vs. 
Poisson  Acquisition,  Model  2 

Number  of  Blue  Direct  Fire  Munitions  used  in  8  days 
Poisson  Acquisition,  Model  2  versus  Classical  Survival,  Model  1 


0  1000  2000  3000  4000  5000  6000  7000  8000  9000  10000 

Number  of  munitions  used  classical  survival  model 
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C.2  Red  Survivors  in  8  Days:  Classical  Survival  Model  1  vs.  Poisson 
Acquisition,  Model  2 

Number  of  Surviving  Red  Platforms  after  8  days 
Poisson  Acquisition,  Model  2  versus  Classical  survival,  Model  1 
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Number  of  munitions  used  classical  survival, 
model  1  with  re-allocation 


D.l  Munitions  Used  by  Blue  in  8  Days:  Classical  Survival,  Model  1,  with 
Proportional  Blue  Reallocation  vs.  Classical  Survival,  Model  1,  No  Blue 
Reallocation 

Blue  Direct  Fire  Munitions  used  during  8  days 
Classical  Survival,  Model  1  versus  Classical  Survival,  Model  1  with  re-allocation 
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Number  of  munitions  used  Poisson  acquisition 
model  with  re-allocation 


D.2  Munitions  Used  by  Blue  in  8  Days:  Poisson  Acquisition,  Model  2  with 
Proportional  Reallocation  vs.  Poisson  Acquisition,  Model  2,  No 
Reallocation 

Blue  Direct  Fire  Munitions  used  in  8  Days 
Poisson  Acquisition  Model  versus  Poisson  Acquisition  Model  with  Re-allocation 


Number  of  munitions  used  Poisson  acquisition  model 
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Number  of  surviving  Red  platforms  for  classical  survival,  model 
1  with  re-allocation 


E.l  Red  Survival  in  8  Days:  Classical  Survival,  Model  1,  with  Proportional 
Blue  Reallocation  vs.  Classical  Survival,  Model  1,  No  Blue  Reallocation 

Number  of  Surviving  Red  platforms  after  8  days 
Classical  Survival,  Model  1  versus  Classical  survival,  Model  1  with  Re-allocation 


0  500  1000  1500  2000  2500 

Number  of  surviving  Red  platforms  for  classical  survival,  model  1 
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Number  of  surviving  Red  platforms  Poisson 
acquisition  model  with  re-allocation 


E.2  Red  Survival  in  8  Days:  Poisson  Acquisition,  Model  2,  with  Proportional 
Blue  Reallocation  vs.  Poisson  Acquisition,  Model  2,  No  Blue  Reallocation 

Number  of  Surviving  Red  Platforms  after  8  days  of  combat 
Poisson  Acquisition  Model,  Model  2  versus  Poisson  Acquisition,  Model  2  with  Re¬ 
allocation 


Number  of  surviving  Red  platforms  for  Poisson  acquisition  model 
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Number  of  surviving  Red  platforms  for  ATCAL- 
like  Model 


F.l  Red  Survival  in  8  Days:  Classical  Survival,  Model  1,  with  ATCAL-like 
Time-Varying  Acquisition  Rate  vs.  Classical  Survival,  Model  1,  No 
Acquisition  Rate  Change 

Number  of  Surviving  Red  Platforms  after  8  days  of  combat 
ATCAL-Like  Model  versus  Classical  Survival  Model  1 
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Number  of  munitions  used  in  ATCAL-like  model 


F.2  Munitions  Used  by  Blue  in  8  Days:  Classical  Survival,  Model  1,  with 
ATCAL-like  Time-Varying  Acquisition  Rate  vs.  Classical  Survival, 
Model  1,  No  Acquisition  Rate  Change 

Blue  Munitions  Used  in  8  Days 
Classical  Survival  Model  1  versus  ATCAL-  Like  Model 
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Appendix  A 

Stochastic  Discrete-Time  Models 

The  previous  pages  of  this  report  describe  a  deterministic  or  “fluid  approximation” 
model,  albeit  between  multi-type  forces.  Here  we  show  how  a  related  stochastic  model 
can  be  formulated.  Further  elaborations  are  possible,  and  will  be  presented  subsequently. 
The  following  stochastic  state  variables  describe  the  situation. 

R{t)  =  Random  variable  denoting  number  of  Reds  at  beginning  (or  selected 
moment)  during  period  t  =  multiple  of  h  >  0:  t  =  0,  h,  2h, . . . ,  1 3h, . . . 

(A.l) 

B(t)  =  Random  number  of  Blues  corresponding. 

Treat  the  co-evolution  of  (R(r),  B(t),  t  =  0,  h,  2h, ... }  as  a  Markov  chain. 

R(t  +  h)  =  R$(t ) 

(A.2) 

B(t  +  h)  -  Bs(t ) 

where  Rs(t),  Bs(t )  are  the  random  numbers  of  survivors  of  period  t  that  are  available  for 
conflict  in  the  subsequent  period  of  duration  (h);  the  state  of  the  system  at  the  beginning 
of  period  t  is  (R(t),  B(t)),  and  it  evolves  to  (Rs(t),  Bs(t ))  at  the  end  of  that  period,  which 
defines  the  state  at  the  beginning  of  period  t  +  h,  and  so  on. 

Poisson  Acquisition  and  Fire 

Consider  the  following  ways  of  assigning  Blues  to  search  for  and  engage  Reds. 

(A)  Sectorization:  Direct  Fire 

Consider  Blue-force  management  first.  At  the  beginning  of  period  t  assign  B(t)/R(t) 
distinct  Blue  shooters,  a  group  of  Blues,  to  each  Red;  e.g.  if  B(t)  =  200  and  R(t)  =  100 
then  the  first  two  Blues  are  assigned  to  the  first  Red,  the  second  two  Blues  are  assigned  to 
the  second  Red,  etc.  If  fractions  arise  treat  them  as  fractions  of  the  time  interval  of 
duration  h. 

Let  Red  forces  engaging  in  sectorized  direct  (aimed)  fire  be  assigned  similarly; 
R(t)/B(t)  distinct  Reds  per  individual  Blue.  Note  that  it  is  not  necessary  that  each  side  be 
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in  the  same  firing  mode.  The  objective  of  “direct  fire”  is  to  assign  an  equal  number  of 
shooters  to  each  target,  and  not  to  allow  random  (e.g.  equally  likely  and  independent) 
target  picking.  This  is  analogous  to  the  “aimed  or  direct  fire”  protocol  that  leads  to  the 
Lanchester  Square  Law. 

Model  the  number  of  contacts/encounters  by  a  Blue  group  with  “its”  Red  during  ( t , 
t  +  h)  as  a  Poisson  process  of  rate  ABR(B(t)/R(t)),  where  the  encounter  rate  parameter  XB 
can  actually  be  a  function  of  the  period,  t,  and  other  variables;  some  can  be  decision  or 
control  variables.  At  each  encounter  with  “its”  Red  let  the  probability  of  kill  be  kbr, 
independently  of  all  other  events.  Therefore  the  kill  process  is  effectively  a  marked 
terminating  Poisson  process ,  and  the  probability  that  a  particular  Red  avoids  being 
killed  -  survives  -  period  t  is  We  make  the  conventjon  that  killed  Reds  are 

removed  at  the  end  of  the  period  (so  they  are  available  to  search  and  shoot  during  the 
period. 

The  same  modeling  assumption  is  adopted  to  describe  a  single  Blue  survivorship 
probability:  *-*"(*(')/*« )*«*_ 

The  above  allows  one  to  write  down  the  one-step  transition  probabilities  of  a  Markov 
chain  {R(t),  B(t),  t  =  0,  h,  2h, ... }  that  is  appropriate  when  the  two  opponents  both  are  in 
direct-fire  mode. 


?{R(r  +  h)  =  r,  B(t  +  h)  =  A|jR(/),  B(f )}  = 


(*?’} 


x 


B(t 


R(i))k  BRh 


9-^RB{R{t)j  B(t)))C  RBh 


\  _  e-*BR{B{,)lR{t))KBRh  r  Q<r< 


R(t) 


(A3) 


l _ e-M*M/*('))^|s(')  b  o<b< 


B(t) 


It  would  be  possible  to  refine  the  above  to  take  account  of  the  order  in  which  Blues 
and  Reds  are  killed  within  ( t ,  t  +  h ),  but  at  the  cost  of  greater  complexity  than  is 
consistent  with  the  present  discrete-time  model  data,  e.g.  from  COSAGE.  We  make  the 
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artificial  convention  that  shot  effects  (kills)  are  resolved  at  the  end  of  an  interval,  so  Reds 
and  Blues  alive  until  the  end  of  period  t  are  alive  and  able  to  shoot  until  (t  +  h)-,  i.e.  the 
instant  before  time  period  t  +  h  starts.  This  convention  could  be  changed,  e.g.  to  have  the 
settle-up  time  be  the  middle  of  an  interval,  i.e.  at  t  +  hi 2,  but  this  seems  unnecessary  and 
is  not  done  here. 

Transition  to  “Continuous  Time”:  Relation  to  Lanchester  “Square  Law” 

It  is  seen  from  the  form  of  each  component  of  the  Markov  transition  probability  (A.3) 
that  the  conditional  means  are 

(')>  5(f)]  =  (A  4) 

and 


(A.6) 


E[B(t  +  h)\R{t),B{t)]  =  £(r)e-M*(')/*(0)*w>  _  (A  5) 

If  h  -*  0,  so  time  intervals  are  short,  then  the  two-term  Taylor  expansion  provides 
E[R{t  +  h)\R(t),  5(f)]  =  R(t)[l  -  XBR  (B(t)/R(t))KBRh  +  0(h 2 )] 
which  leads  to  the  differential  equation 

and  (formally)  removing  conditions  yields 


d_ 

dt 


(A.7) 


(A. 8) 


which  are  precisely  the  form  of  the  Classical  Lanchester  Square  Law  equations.  The 
boundary  conditions  £[/?(/)]  >  0  and  E[B(t)]  >  0  must  be  imposed  for  correct  solution. 

Consequently,  we  argue  that  the  deterministic  fluid  equations  given  earlier  are  a 
consequence  of  this  stochastic  model,  and  an  approximation  to  aspects  of  it.  Furthermore, 
variability  can  be  assessed  straightforwardly: 
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and  the  entirely  analogous  equation  for  Var^B(t  +  /z)|  R(r),  B(r)] . 

Firing  Options  and  Immediate  (Imperfect)  BDA 

The  model  proposed  can  be  easily  implemented  to  study  several  firing  options: 
repeated,  information  affected  shots  at  the  same  target;  they  can  do  so  under  the  guidance 
of  local  BDA  conducted  soon  after  a  shot.  Furthermore,  the  BDA  can  be  rendered  as 
realistically  error-afflicted.  Incorporation  of  the  various  random  effects  (shot  hits  and 
misses,  BDA  assessments,  shot  repetitions,  repeated  BDA,  etc.). 

For  the  present  we  model  BDA  effects  at  the  unit  engagement/shot  level.  The 
information  presumed  used  to  guide  subsequent  shots  is  obtained  at  the  time  of  shooting 
(just  following  a  shot);  it  may  be  provided  by  the  shooter,  e.g.  from  the  actual  shooting 
ground  unit,  from  a  forward  observer,  or  possibly  from  an  air  observer  (UAV  or  helo)  — 
or  combination  thereof.  Note  that  for  what  is  done  here  the  observers  are  assumed  present 
and  connected  to  the  shooter.  This  may  not  be  so,  but  the  effect  of  loss  can  be  modeled  in 
various  ways  (not  included  here). 

Here  are  some  options. 

(a)  Every  engagement  involves  a  single  shot,  e.g.  by  Blue  or  Red.  In  this  familiar 
case  we  put  kbr  =  kBr{  1 )  in  basic  attrition  equations  (A.3);  here  *&«(1),  or  /&b(1),  is  a 
single-shot  kill  probability.  It  should  be  regarded  as  conditional  on  as  many  influential 
variables  as  are  accessible. 

(b)  Every  engagement  requires  or  is  ordered  to  fire  (by  doctrine)  a  single  salvo  of  sB 
shots.  Then,  simply,  assuming  independence,  the  effective  engagement  kill  probability  is 
kbr  =  1  -  (1  -  Kb/?(1))  .  Or  there  may  be  a  pattern  fired;  if  symmetrically  placed,  this 
might  be  equivalent  to  kbr  =  1  -  (1  -f  kbr{  1))  (1  -  f2icBR(  1))  for  sB  =  2,  where  f  (>  0)  and 
fi  0)  are  selected  to  “automatically”  adjust  to  possible  target  movement.  Again, 
substitute  into  (A.3)  or  the  equivalent  and  solve  the  equations. 
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(c)  Suppose  engagements  are  quick-time  sequences  of  Shoot,  Look,  Shoot  (SLS) 
with  some  termination  rule  (e.g.  second  shot  is  last).  The  “shot”  is  then  a  nearly 
instantaneous  random  sequence  of  individual  shots  punctuated  by  brief  effect 
assessments,  which  are,  unfortunately,  potentially  incorrect: 

CBR(k)  =  P{Blue  shot  estimated  to  have  killed  Red  target  I  Blue  shot  killed  Red 
target,  or  Red  target  dead;  k  means  actually  killed} 

cbr(™)  =  PfBlue  shot  estimated  to  have  killed  Red  target  I  Blue  shot  missed  Red 
target,  or  Red  target  alive;  m  means  actually  missed/alive} 

Other  needed  probabilities  can  be  obtained  by  complementation;  let 
cBR{h)  =  1  -cBR(h),  cBR(m)=l-cBR(m) .  Now  the  resultant  kill  probability  of  (S,L,S)  is 

kBr=kbr(  l)cB/?  (/z)  +  (l  -  kbr  ( l))cB*  (m)  kbr  ( l) ; 

the  second  shot  is  here  assumed  to  have  the  same  kill  probability;  it  need  not.  The 
expected  number  of  shots  of  Blue  vs.  Red  per  engagement  (kill,  or  not)  is,  under  the 
doctrine  (S,L,S), 

e[”br  }  =  1  •  kBr  ( l)cB*  (h)  +  2(l  -  kbr  ( l)) . 

Another  measure  of  effectiveness  is  the  expected  number  of  extra  shots  taken  after  a  kill, 
(if  by  Blue,  eBR)  during/on  a  given  engagement.  This  is,  of  course,  the  consequence  of 
imperfect  BDA.  In  the  S,L,S  case  this  is 

E\_eBR\  =  CBR{h) 

where  this  is  just  the  probability  that  the  second  shot  is  unnecessary. 

(B)  Indirect/Area  Fire 

Suppose  Blue  conducts  Indirect  or  Area  Fire  at  Red  units  in  a  region.  Model  as 
follows:  a  single  Blue  shooter’s  “shot”  (possibly  a  pattern  or  volley)  makes  a  kill  on  each 
Red  unit  independently  with  probability  T]Br ,  so  with  probability  JjBR=l-7jBR  a 

particular  Red  unit  survives.  It  follows  that  if  each  Blue  shoots  according  to  a  Poisson 
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process  with  rate  %br  that  the  particular  Red  unit  survives  that  one  Blue’s  shot  with 
probability  e~^BRh7?BR ,  and  survives  all  of  Blue’s  area  fire  during  (t,  t  +  h )  with  probability 
e  $BRhTlBRB(l)  _  Since  (by  assumption)  all  Reds  are  equivalently  available,  the  probability 
distribution  of  surviving  Reds  is  binomial 


p{R(/  +  /i)  =  r|B(r),R(r)}  =  ^  _*-*«**«*(') j*(,)  ' 


Under  the  previous  assumptions,  e.g.  of  (A),  a  similar  expression  holds  for  B(t  +  h ).  The 
conditional  one-step  transition  probabilities  are  independent,  so  we  arrive  at  system 
transition  probability  expressions  quite  analogous  to  that  for  (A.3)  above.  The  one-step 
transition  probabilities  appropriate  if  each  side  is  exchanging  area  fire  is 

/>{*(? +h)  =  r,  B(t  +  h)  =  fc|R(f),  tf(r)}  = 


Note:  It  is  perfectly  possible  to  have  one  side.  Blue,  say,  attack  Red  using  area  fire,  while 
Red  attacks  Blue  using  direct  fire,  or  to  start  with  all  indirect  fire,  and  transition  to  direct 
fire  as  time  increases  and  range  decreases.  Combinations  of  various  sorts  can  be  modeled, 
and  the  results  analyzed. 
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Appendix  B 

Variability/Uncertainty  Assessment 

In  actual  COSAGE  evolutions  there  are  several  replications  (e.g.  r=  16)  made,  and 
the  results  eventually  averaged  to  get  what  we  have  called  the  data  values  k(jR,jB; 2) 

mdk(jB,jR;2).  These  can  be  directly  used  to  point-estimate  the  probability  of  survival 
k (x,  y;  1 ),  and  also  R(jR:  t),  B(jB;  t),  t  =  0,  h,2h, llh, ... 

In  the  most  accessible  documentation  the  kill  data  is  summarized  by  mean  and 
(estimated)  variance: 

k(x,  y  ;2)  =  -  £  *. ;  ( x,  y  ;2)  =  k  (x,  y  ;2)  (B .  1  ,a) 

^  /=! 

var[fc(x,  y-2)]  =  -L £  [k,  (x,  y;2)  -  k(x,  y;2)f .  (B.  1  ,b) 

r  1  i=i 

Similar  data  are  available  for  the  number  of  weapons  fired. 

A  simple  computational  way  of  checking  for  the  stability/sampling  variability  of  the 

results  used  for  predicting  (mean)  survivors  on  each  side  is  to  (i)  re-sample  each  k-value 
from  a  normal  or  Gaussian  distribution  N(k ,  var[k(x,  y;2)]) ,  where  r  is  the  number  of 

independent  and  identical  (by  assumption)  replications  (e.g.  16),  and  use  the  result(s)  to 
estimate  tc(x,y;2h )  from  (4.10)  and  (4.18);  then  (ii)  apply  this  outcome  as  survival 

parameters  in  equations  (4.10)  and  (3.3  a&b);  repeat  this  independently  ~  50  times  at  t  = 
4,  8,  10,  12  and  examine  the  results  (compute  R(jR-,mh),  B(jB;lh),  m,l=  0,  1, 

2,  ...,  31, ...,  and  the  respective  variances  and  standard  deviations).  This  approach  is  a 
way  of  quantifying  the  uncertainty  in  a  mean  or  deterministic  approximation.  It  may  be  of 

interest  to  use  as  a  “certainty  equivalent”  as  part  of  the  process  a  lower  confidence  level 
(roughly  R(jR',t;q')-2-Jvar /?  ,  or,  possibly  better,  a  value  obtained  by  bootstrapping. 


48 


Appendix  C 

Models  for  Direct  and  Indirect  Fire 

For  simplicity  assume  there  is  one  type  of  Red  target  and  one  type  of  Blue  platform 
firing  one  type  of  Blue  weapon.  Let  5(0)  be  the  initial  number  of  Red  targets  and  5(0)  be 
the  initial  number  of  Blue  platforms.  COSAGE  summary  data  contain  the  averages  of 
various  measures  over  16  replications  of  simulations  for  2  days  of  combat  of  Red  and 
Blue  combatants  in  various  postures. 

Aimed  (or  Direct)  Fire  Weapons 

For  aimed  (direct  fire)  fire  weapons,  COSAGE  summary  data  contain  the  average 
number  of  kills  and  the  average  number  of  shots  fired  over  the  r-  16  replications.  Let 
5(0),  (respectively  5  (0) ),  be  the  total,  (respectively  average),  number  of  shots  over  the  16 
replications.  5(0)  =  165(0).  Let  AT(0),  (respectively  K  (0)),  be  the  total,  (respectively 
average),  number  of  kills  over  the  16  replications.  Each  aimed  fire  weapon  is  shot  at  one 
Red  target.  Let  A(f)  be  equal  to  1  if  the  jth  shot  kills  the  target  and  0  otherwise. 

so 0) 

w)  =  2>0) 

E[K(0)]=S(0)Kw(d) 

Thus  an  estimate  of  the  probability  of  kill  for  an  aimed  fire  weapon  is 

5(  0)  5(0) 


Now  to  avoid  having  to  estimate  a  kill  probability  as  0  owing  to  small  sample  luck  we 


replace  k (d)  =  — — 
5(0) 


g(0) 

5(0) 


by  a  Bayes  estimate  that  assume  kills  are  binomial  with  5 


trials  and  applies 
5(0) +  2 


a  uniform 
g(0)  +  Q/r) 
5(0)  +  (2/ r) 


prior;  the  result  is  the  estimate  of  kill 
where  r  is  the  number  of  replications,  e.g.  16. 
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To  estimate  the  probability  a  particular  Red  survives  all  5  shots,  we  assume  that  the 

number  of  shots  fired  at  a  particular  Red  over  the  r=  16  replications  is  -  ^  =  —  °^ 

rR  R 

Thus,  the  probability  a  particular  Red  survives  all  the  shots  is 

nS(0)« 


Ud)= 


x  5(0) +  (1  /r) 
S(0)  +  (2/  r)_ 


The  rate  of  shooting  per  Blue  over  the  48-hour  period  of  COSAGE. 

m  ’  B(0) 

An  equation  to  compute  the  average  number  of  Reds  that  survive  aimed  (direct)  fire 
during  the  time  interval  (t,  t  +  h]  is 


R(t  +  h)=R(t)icR(d)  = 


=4-1 


g(0)  +  (l/r)' 
5(0)  +  (2 /r) 


p{d  )hB(t)/ R(t) 


An  estimate  of  the  average  number  of  aimed  fire  weapons  shot  by  one  Blue  platform 
against  Red  targets  is  5  (0)/  5(0) .  Assuming  an  infinite  supply  of  Red  targets,  an  estimate 

of  the  expected  rate  of  kill  by  one  Blue  platform  against  Red  targets  is 


Indirect  or  Area  Fire  Weapons 

COSAGE  summary  data  record  the  average  number  of  Red  targets  that  are  affected  by 
the  indirect  fire  shots.  Let  Nj{ 0)  (respectively  N j  (0))  be  the  total  (respectively  average) 

number  of  Red  targets  affected  by  the  5(0)  (respectively  5(0))  indirect  (area  fire) 
(respectively  average)  fire  shots.  Nf(  0)  =  rNf  (0)  =  16A^(0)  Let  /,(/;  a;f)  be  equal  to  1  if 

the  /th  indirect  weapon  shot  affects  the  fh  Red  target  and  be  equal  to  0  otherwise. 

>=1  y=i 
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Let  Xt  be  equal  to  1  if  the  j*  Red  target  affected  by  an  indirect  weapon  shot  is  killed  by 
the  weapon  and  0  otherwise.  The  number  of  Red  targets  killed  by  indirect  (area)  weapon 
shots  is 


Nf(0)  5(0)  R(0) 

w=X^  =  II'iw) 

/=!  /=1  j=\ 

where  /,(/;  a;  k)  is  equal  to  1  if  the  ith  indirect  (area)  weapon  shot  kills  the  fh  Red  target 
and  is  0  otherwise. 


Nf(0)  s  R 

E[jf(0)]=  E[  £  X,]  =  £[££/, =  smm.(«) 

i= 1  ;=1  j- 1 

Thus  an  estimate  of  the  probability  an  indirect  weapon  shot  kills  a  Red  target  is 

it  (3)_  ^(°)  _  Jm 

^  S(0)R(0)  S  (0)R(0) 

where  ^T(O)  is  the  average  number  of  Red  targets  killed  over  the  r  COSAGE  replications. 

To  avoid  having  to  estimate  a  kill  probability  as  0  owing  to  small  sample  luck  we 

.  .  /  .  AT(0)  K( 0)  ^  . 

replace  K  \a)  = - =  ^= -  by  a  Bayes  estimate  that  assumes  kills  are 

S(0)R(0)  S  (0)R(0) 


binomial  with  S(0)R(0)  trials  and  applies  a  uniform  prior;  the  result  is 


i  (_Tt)-  _&(S»  +  (Ur) 

’  ^  S(0)R(0)  +  2  S(0)R(0)  +  {2/r) 

Estimate  the  number  of  shots  taken  over  a  single  COSAGE  replication  as  S( 0) .  The 
probability  a  particular  Red  survives  all  S  (0)  indirect  shots  is  estimated  as 


_g(0)  +  (l/r) 
S(0)R(0)  +  (2/r) 


The  average  number  of  indirect  fire  shots  fired  per  Blue  over  the  initial  48-hour 
period  of  COSAGE  is 
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An  equation  to  compute  the  average  number  of  Reds  surviving  all  the  indirect  fire 


shots  in  (t,  t  +  h]  is 


R(t  +  h)  =  R(t)K(a)  =  R(t) 

'  _  S  (0)R(0)  +  (2  /  r) 


p{a)hB(t) 


The  expected  number  of  Reds  killed  by  indirect  fire  shots  during  ( t ,  t  +  h]  is 


The  average  number  of  indirect  weapons  fired  by  each  Blue  at  Red  targets  is 
5(0)/  5(0) .  Assuming  an  infinite  supply  of  Red  targets,  an  estimated  average  rate  at 

which  one  Blue  platform  using  area  fire  weapons  kills  Red  targets  is  during  the  first  48 
hour  period  modeled  by  COSAGE  is 

sia1o) = MJmifcL = ffa)-if(0)+(|/r> 

B(0)  S  (0)R(0)  +  (2/r)  V  S(0)R(0)  +  {2/ r) 

The  estimated  average  rate  at  which  one  Blue  platform  using  area  fire  weapons  kills  Red 
targets  during  ( t ,  t  +  h]  is 


Example  1: 


S(a;t)=  p(a)h 


_  AT(0)  +  (l/ r ) 
S(0)R(0)  +  (2/r) 


UH155Z  firing  area  fire  munition  M483A1  at  platform  R1NTP 

Given:  Initial  number  of  RINTP  =  1 1068 
Initial  number  of  UH155Z  =  144 

Average  total  number  of  indirect  shots  SBR  {jB,wB)  =  263 1 8 
Observed  number  of  kills  KBR  ( jB  ,wB,jR )  =  509.56 
Number  of  Replications  =16 


Bayes  probability  of  kill  for  a  single  shot 


„  /  ,x  K  + 1  K  +  (l/r) 
rja;b)  =  — —  =  \  =  • 


509.56 +  - 


16 


SR  +  2  SR  +  {2/r)  (26318*1 1068)  + (2/16) 

The  probability  a  Target  survives  all  shots  is 


=  1.74955 £-06 


?*(«)  = 


\  K  +  (l/r)' 
SR  +  {2/r) 


=  [l  - 1 .749552s  -  06f318  =  .954999 


The  estimated  average  number  killed  in  48  hours  is 

(#  targets)  x  [l  -  qR (a)]  =  1 1068  x  (1-.954999)  =  498.0682 

which  is  close  to  the  COSAGE  number  of  509.56. 

Example  2: 

UH155Z  firing  area  fire  munition  M483A1  at  platform  ROPR7 

Given:  Initial  number  of  ROPR7  =  7560 
Initial  number  of  UH155Z  =144 

Average  total  number  of  indirect  shots  SBR  ( jB ,  wB )  =26318 
Observed  number  of  kills  KBR  ( jB ,  wB ,  jR )  =41 6.62 
Number  of  Replications  =16 

Bayes  probability  of  kill  for  a  single  shot 


£w{a-,b)  = 


K  + 1  _  2sT  +  (l/r) 


416.62  +  - 


16 


SR  +  2  SR  +  {2/r)  (26318*7560)+ (2/16) 


2.09426is  -  06 


The  probability  a  Target  survives  all  shots  is 


\  K+{l/r) 

SR  +  {2/r) 


[l  -  2.09426E  -  06]26318  =  .946375 


The  estimated  average  number  killed  in  48  hours  is 


(#  targets) x [l - qR(a)]  =  1 1068x  (1-.946375)  =  405.4079 
which  is  close  to  the  COSAGE  number  of  416.62. 
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Appendix  D 

ATCAL-Like  Firing  Rate  Model 

Versions  of  ATCAL  (abbreviation  of  An  Attrition  Model  Using  Calibrated 
Parameters)  Modeling  in  DISC-O-TIC  Style 

About  twenty  years  ago,  in  the  early  1980s,  the  Army,  through  efforts  at  the  then 
Combat  Analysis  Agency  (now  Center  for  Army  Analysis)  or  CAA,  developed  An 
Attrition  Model  Using  Calibrated  Parameters ,  or  ATCAL:  “a  new  method  for  calibrating 
a  set  of  attrition  equations  to  the  results  of  sample  high-resolution  simulations.”  The 
algorithm  developed  is  a  set  of  attrition  equations  for  “point”  (we  call  it  direct )  and  “area” 
(we  call  it  indirect )  fire;  their  algorithm  (equations)  computes  losses  by  cause  (round,  or 
ammunition  type),  using  high-resolution  simulation  data,  as  we  have  done  earlier  in  this 
paper;  thus  ATCAL  parameters  are  estimated  from  high-resolution  data,  as  are  ours. 
ATCAL  does  recognize  target  “availability”  e.g.  visibility,  in  a  simplified  probabilistic 
manner,  and  also  target  priorities:  such  priorities  are  governed  by  intrinsic  importance  to  a 
shooter,  but  also  are  higher  for  those  more  easy  to  kill  than  other  available  targets.  For  the 
present,  DISC-O-TIC  does  not  go  to  such  explicit  lengths  because  of  limits  in  COSAGE 
data  available.  However,  ATCAL  does  not  attempt  to  model  the  time-step-evolution  or 
dynamics  of  combat,  whereas  DISC-O-TIC  does,  as  seen  above.  ATCAL  seems  to 
assume,  rather  specifically,  that  the  forces  decrease  exponentially  with  the  duration  of  a 
battle,  which  is  an  effect  that  occurs  if  one  reasonably  assumes  that  fire  allocation  to  a 
target  type  should  decrease  as  the  number  of  such  targets  is  decreased. 

Here  is  a  way  of  adjusting  firing  (actually  acquisition  rate)  in  DISC-O-TIC  to  achieve 
the  above  effect. 

Replace  the  acquisition  rate  per  shooter  in  Classical  Survival  Model  1,  nsQ(*,js,ws 5  jo) 
by  a  time-dependent  rate,  n$Q(*,  js ,w$;  jQ,t) ;  here  S  denotes  Shooter  type,  and  Q  is 
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quarry/target  type;  if  S  is  Red  then  the  Q  is  Blue,  and  vice  versa.  Total  shooting  rate  per 
quarry  is,  after  replacement 

nsQ  (•,  js ,  ;  Jq  ,t)S(js,wS,  t)l  Q(  jQ ,  t )  (D.  1 ) 

total  killing  rate  per  quarry  is 


n*sQ (•,  js , ; ]q , t) (ln(xr5e ))s {js,ws,t)/Q(jQ,t) .  (D.2) 

Now  suppose  the  shooting  rate  is  actually  modified  by  the  factor  in  square  brackets: 

QUd’t) 


nsQ{*,js,ws;jQ,t)  = 


S(js,ws,t ) 


/sqUs’WsJq) 


nSQ(*,js,Ws‘jQ ) 

I  V - - - ' 

! shooting  rate  deduced  from 
correction  term  that  reduces  (0*2 h]  COSAGE  data  and 

total  Shooting  rate  if  Direct  Fire  Model  1 

Quarry/Shooter  ratio  decreases 


(D.3) 


-  nSQ (•,  js , W5 ;  ]q )fsQ (js , ws ;  jQ ) ~JrrQ  ^ 


S(js,ws,t)’ 

fsQ(js,  jo)  is  here  a  constant  to  be  determined. 

This  particular  form  of  the  Shooter-Quarry  rate,  if  applied  across  all  shooter  options 
leads  to  geometric  decline  of  the  Quarry  forces  in  time  although  not  necessarily  at  the 
same  rate  by  type.  Prioritization  by  type  can  be  achieved  by  adjustment  of  the  tuning 
constant  JSq:  it  can  be  replaced  by  fsQ(js,  jQ). 

Given  the  shooting  rate  per  period  we  get  for  Model  1  and  the  ATCAL-like  rule. 


Q(jQ,t  +  h)  =  Q(jQ,t)exp 


X  X nsQ  (*’  J‘s » ' ;  Jq  ) ln  Ksq(j's  ,  w5 ;  jQ  )fsQ{js ,  ws ;  jQ ) 


JS  W5 


•(D.4) 


=  QUq^CqUq) 


so 


Q(Jq>  0  -  QUq^)[cq(Jq)] 


(D.5) 
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To  compute  the  expenditure  of  vvs-type  weapons/munitions  up  to  the  end  of  time  period  t 
it  is  necessary  to  sum  expenditures  during  each  period:  if  M(wQ ;  t)  is  the  expenditure  of 
ws  through  period  t. 


M(wQ ; ' 0  =  2 2 £  ■ n*SQ (•, js , Ws ;  Jq , X ) ^(-’X\  Q{jQ,t) 
is  jQ  Q{JQ,x) 


hsq{*’  js ,  ws ;  jo)fsQ{js » ws ;  jo)Q{jQ ,  x) 

X=i  js  }Q 


(D.6) 


_  / 

=  /L, X nSQ (•’ js > Ws ; ]Q )f SQ ( js , ws ; jQ )  £ Q(jQ,x ) 

.  rs  )Q  JC=1 

The  function  fsaijs,  ws'Jq)  is  a  potentially  arbitrary  control  function.  It  follows  from 
(D.5)  that  one  way  to  determine  its  value  is  to  use  the  initial  COSAGE  data;  let  the 
COSAGE  Shooter-weapon-specific  (js ,  ws)  number  of  survivors  of  Quarry  type  Jq  be 
q(js,ws;jQ,  2 h). 

Then,  replacing  fSQ  by  its  estimate  JSq  from  available  COSAGE  data, 

q(js,ws;jQ,2h)  =  q(js ,  w5 ;  jQ  ,0)  exp[n(*,  js ,  ;  jQ  )  In  KSQfSQ  ( js ,  w5 ;  jQ)]  .  (D.7) 

Now  solve  for  the  initial-data-driven  value  f$Q : 


fsQ(js,ws;jQ)  =  -r-: - :  *_  (D8) 

'"[•>Js,wsljQ)lnKsQ{js,ws;jQ)  [q(js,ws;jQ, 0)1 


In  the  special  case  when 


(a)  q(js,ws;jQ,2h)  =  q(js,ws;jQ,  0) ,  then  we  put  fsQ{js,ws;jQ)  =  S(js,  0); 


(^)  q(js,ws',  jQ ,2h)  —  0 ,  then  we  put  /sq ( js ,ws~,  jo)  —  S[ js ,0)  for  the  initial  time 
period  and  then  0  for  the  other  time  periods;  and  the  expenditure  of  munitions 


M{js,ws-JQ,t) 


\nsQ  (•»  js ,  ws ;  jQ  )S  ( js  ,0) 

jo 


for  t-2h 
for  /  >  3  h 
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Discussion.  Our  basic  survival  model  (D.4)  and  (D.5)  is  a  consequence  of  the  Shooter 
being  able  to  perceive  his  own  and  the  quarry  force  types  and  sizes  perfectly  at  each  time 
which  is  manifestly  highly  optimistic.  Sensitivity  to  this  assumption  can  be  studied  if  only 
by  simulation.  This  is  left  for  future  work. 

No  particular  reason  has  been  given  for  adjusting  acquisition/shooting  rate  to 
achieve  (nearly  precise)  geometric/exponential  decrease,  presumably  but  not  necessarily 
on  both  sides.  This  choice  does  drive  the  quarry  force  size  down  while  limiting  munition 
expenditures,  but  there  are  other  variations  possible  in  the  acquisition/firing  rates  worth 
exploration;  one  or  more  of  these  might  provide  superior  results  from  some  standpoint; 
they  should  be  investigated  and  clarified  in  later  work. 

If  more  detailed  output  from  COSAGE  (or  other  such  high-resolution  models) 
were  available,  it  could  be  possible  to  add  additional  detail  to  our  meta-model, 
DISCOTIC,  that  would  provide  greater  insight  into  choice  of  suitable 
acquisition/shooting  rates.  Also,  the  cost  effectiveness  of  surveillance,  e.g.  by  helicopter, 
could  be  studied  by  meta-modeling.  Such  studies  could  effectively  guide  future  high- 
resolution  model  (e.g.  COSAGE,  but  not  exclusively)  runs,  and  provide  the  present 
models  with  greater  modem  detail. 
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Appendix  E 


Weapon  reallocation 

Weapon  platforms  will  fire  the  same  type  of  munitions  at  several  types  of  targets.  As 
target  platforms  are  attrited  during  combat,  weapon  platforms  will  re-allocate  their  fires 
among  the  surviving  platforms.  One  perhaps  overly  simple  re-allocation  algorithm  is  to 
reallocate  based  on  the  original  proportion  of  one  kind  of  weapon  fired  at  each  type  of 
target;  Washburn  (2000).  In  particular,  the  number  of  weapons  of  type  wB  fired  by 
platforms  of  type  jB  at  targets  of  type  jR  during  the  time  interval  (t,  t+h]  under  this  re¬ 
allocation  algorithm  is 


NBR{jB,wB,jR,t  +h)~  NBR[jB,wB,jR,t) 


1  ^BR{jB^B,jR2h)l[R{jB,t)  >  l) 


where  l(R(jB,t}>  l)— 1  if  R(jBj} >  1  and  0  otherwise.  This  reallocation  algorithm  has 


the  effect  of  focusing  the  firing  of  weapons  of  this  type  on  the  remaining  target  platforms. 
It  ignores  possible  delays  in  target  acquisition  as  the  number  of  target  platforms  dwindles. 


58 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center . 2 

8725  John  J.  Kingman  Rd.,  Ste  0944 

FT.  Belvoir,  VA  22060-6218 

2.  Dudley  Knox  Library,  Code  013 . 2 

Naval  Postgraduate  School 

Monterey,  CA  93943-5100 

3.  Research  Office,  Code  09 . 1 

Naval  Postgraduate  School 

Monterey,  CA  93943-5138 

4.  Prof.  Donald  P.  Gaver  (Code  OR/Gv) . 1 

Dept  of  Operations  Research 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

5.  Prof.  Patricia  A.  Jacobs  (Code  OR/Jc) . 1 

Dept  of  Operations  Research 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

6.  Dr.  Paul  K.  Davis . 1 

RAND  Corporation 

1700  Main  St. 

Santa  Monica,  CA  90406 

7.  Dr.  Bruce  W.  Fowler . 1 

USAAMCOM 

ATTN:  ANSAM-RD-AS 
Redstone  Arsenal,  AL  35898-5242 

8.  Mr.  Walter  W.  Hollis . 1 

Deputy  Under  Secretary  of  the  Army  (OR) 

ATTN:  SAUS  (OR) 

The  Pentagon,  Room  2E660 
Washington,  DC  2031 0-0 1 02 


59 


9.  Dr.  Moshe  Kress . 

1 

CEMA 

P.O.B.  2250  (TI) 

Haifa  31021 

ISRAEL 

10.  Dr.  Mark  A.  Youngren . 

1 

2929  Grey  St. 

Oakton,  VA  22124 

1 1 .  LTC  J.  Appleget . 

1 

TRADOC-Monterey 

Naval  Postgraduate  School 

166  Bouldry  Rd,  Rm  212 

Monterey,  CA  93943-5213 

12.  Prof.  Kevin  R.  Wood  (Code  OR/Wd) . 

Dept  of  Operations  Research 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

. 1 

13.  Prof.  Tom  Lucas  (Code  OR) . 

Dept  of  Operations  Research 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

. 1 

14.  Mr.  Michael  Bauman . 

1 

Director, 

USA  Training  &  Doctrine  Command  Analysis  Center  (TRAC) 

Fort  Leavenworth,  KS  66027 

15.  Mr.  John  Warren . 

1 

U.S.  Army  Concepts  Analysis  Agency 

8120  Woodmont  Avenue 

Bethesda,  MD  208 1 4 

16.  Mr.  H.  Kent  Pickett . 

. 1 

Director,  Modeling  and  Research  Directorate 

USA  Training  &  Doctrine  Command  Analysis  Center  (TRAC) 

Fort  Leavenworth,  KS  66027 


60 


17.  Mr.  M.  Mahony . 1 

U.S.  Army  Concepts  Analysis  Agency 

8120  Woodmont  Avenue 
Bethesda,  MD  20814 

18.  Mr.  Yee  Kah  Chee . 1 


Senior  Manager  (Operations  Research) 

DSO  National  Laboratories 
20,  Science  Park  Drive 
SINGAPORE 

19.  Mr.  Yeo  SiokKhoon . 1 

Operations  Research 
DSO  National  Laboratories 
20,  Science  Park  Drive 
SINGAPORE 


20.  Prof.  Gordon  Schacher . 1 

Institute  for  Joint  Warfare  Analysis 

Naval  Postgraduate  School 
Monterey,  CA  93943-5101 

21 .  Mr.  Curtis  Blais . 1 

Institute  for  Joint  Warfare  Analysis 

Naval  Postgraduate  School 
Monterey,  CA  93943-5101 

22.  Prof.  William  Kemple . 1 

Institute  for  Joint  Warfare  Analysis 

Naval  Postgraduate  School 
Monterey,  CA  93943-5101 

23.  Prof.  Alex  Callahan..- . 1 

Institute  for  Joint  Warfare  Analysis 

Naval  Postgraduate  School 
Monterey,  CA  93943-5101 

24.  Prof.  Shelley  Gallup . 1 

Institute  for  Joint  Warfare  Analysis 

Naval  Postgraduate  School 
Monterey,  CA  93943-5101 


61 


